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HIGHER-ORDER ELLIPTIC EQUATIONS IN NON-SMOOTH 
DOMAINS: HISTORY AND RECENT RESULTS 

ARIEL BARTON AND SVITLANA MAYBORODA 


Abstract. Recent years have brought significant advances in the theory of 
higher order elliptic equations in non-smooth domains. Sharp pointwise es¬ 
timates on derivatives of polyharmonic functions in arbitrary domains were 
established, followed by the higher order Wiener test. Certain boundary value 
problems for higher order operators with variable non-smooth coefficients were 
addressed, both in divergence form and in composition form, the latter being 
adapted to the context of Lipschitz domains. These developments brought new 
estimates on the fundamental solutions and the Green function, allowing for 
the lack of smoothness of the boundary or of the coefficients of the equation. 
Building on our earlier account of history of the subject in [ 25], this survey 
presents the current state of the art, emphasizing the most recent results and 
emerging open problems. 
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1. Introduction 


The theory of boundary value problems for second order elliptic operators on Lip¬ 
schitz domains is a well-developed subject. It has received a great deal of study in 
the past decades and while some important open questions remain, well-posedness 
of the Dirichlet, Neumann, and regularity problems in L p and other function spaces 
has been extensively studied in the full generality of divergence form operators 
— divAV with bounded measurable coefficients. 

The corresponding theory for elliptic equations of order greater than two is much 
less well developed. Such equations are common in physics and in engineering 
design, with applications ranging from standard models of elasticity 11021 to cutting- 
edge research of Bose-Einstein condensation in graphene and similar materials 1125 . 
They naturally appear in many areas of mathematics too, including conformal 
geometry (Paneitz operator, Q-curvature [31], [32]), free boundary problems [I], 
and non-linear elasticity |134|. m , [ai- 

it was realized very early in the study of higher order equations that most of the 
methods developed for the second order scenario break down. Further investigation 
brought challenging hypotheses and surprising counterexamples, and few general 
positive results. For instance, Hadamard’s 1908 conjecture regarding positivity of 
the biharmonic Green function [57] was actually refuted in 1949 (see [47] . 55! . 
[127] ). and later on the weak maximum principle was proved to fail as well, at least 
in high dimensions [ 100] , [120 ,. Another curious feature is a paradox of passage to 
the limit for solutions under approximation of a smooth domain by polygons m, 


[95]. 


For the sake of concreteness, we will mention that the prototypical example of 
a higher-order elliptic operator, well known from the theory of elasticity, is the 
bilaplacian A 2 = A (A); a more general example is the polyharmonic operator A m , 
m > 2. The biharmonic problem in a domain fl C R” with Dirichlet boundary 
data consists, roughly speaking, of finding a function u such that for given /, g, h, 


( 1 . 1 ) 
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subject to appropriate estimates on u in terms of the data. To make it precise, as 
usual, one needs to properly interpret restriction of solution to the boundary u\ gn 
and its normal derivative d v u\ gn , as well as specify the desired estimates. 

This survey concentrates on three directions in the study of the higher order 
elliptic problems. First, we discuss the fundamental a priori estimates on solutions 
to biharmonic and other higher order differential equations in arbitrary bounded 
domains. For the Laplacian, these properties are described by the maximum prin¬ 
ciple and by the 1924 Wiener criterion. The case of the polyharmonic operator has 
been only settled in 2014-2015 85], ]84], and is one of the main subjects of the 
present review. Then we turn to the known well-posedness results for higher order 
boundary problems on Lipschitz domains with data in L p , still largely restricted 
to the constant coefficient operators and, in particular, to the polyharmonic case. 
Finally, we present some advancements of the past several years in the theory of 
variable coefficient higher order equations. In contrast to the second order opera¬ 
tors, here the discussion splits according to severals forms of underlying operators. 
Let us now outline some details. 

On smooth domains the study of higher order differential equations went hand- 
in-hand with the second order theory; in particular, the weak maximum principle 
was established in 1960 ([6]; see also Il03l . i 1041 1. Roughly speaking, for a solution 
u to the equation Lu = 0 in 11, where L is a differential operator of order 2m and 
f l C K” is smooth, the maximum principle guarantees 

(1.2) max ||9“u|| i oo (n) < C max ||9^u|| Loo(an) , 

with the usual convention that the zeroth-order derivative of u is simply u itself. 
For the Laplacian (m = 1), this formula is a slightly weakened formulation of the 
maximum principle. In striking contrast with the case of harmonic functions, the 
maximum principle for an elliptic operator of order 2 to > 4 may fail, even in a Lips¬ 
chitz domain. To be precise, in general, the derivatives of order (to— 1) of a solution 
to an elliptic equation of order 2 to need not be bounded. However, in the special 
case of three dimensions, was proven for the m-Laplacian (—A) m in domains 
with Lmschitz boundary. ('(1181. 11201: see also m, ma, ra, m\ for related 
work), and, by different methods, in three-dimensional domains diffeomorphic to a 
polyhedron ('[75]. IfOO]). 

Quite recently, in 2014, the boundedness of the (to — l)-st derivatives of a so¬ 
lution to the polyharmonic equation (—A ) m u = 0 was established in arbitrary 
three-dimensional domains [85]. Moreover, the authors derived sharp bounds on 
the /c-th derivatives of solutions in higher dimensions, with k strictly less than to — 1 
when the dimension is bigger than 3. These results were accompanied by pointwise 
estimates on the polyharmonic Green function, also optimal in the class of arbi¬ 
trary domains. Furthermore, introducing the new notion of polyharmonic capacity, 
in [84j the authors established an analogue of the Wiener test. In parallel with 
the celebrated 1924 Wiener criterion for the Laplacian, the higher order Wiener 
test describes necessary and sufficient capacitory conditions on the geometry of 
the domain corresponding to continuity of the derivatives of the solutions. Some 
earlier results were also available for boundedness and continuity of the solutions 
themselves (see, e.g., EZ), m, EH)- 

We shall extensively describe all these developments and their historical context 
in Section [3] 
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Going further, in Sections |4] and [5] we consider boundary value problems in ir¬ 
regular media. Irregularity can manifest itself through lack of smoothness of the 
boundary of the domain and/or lack of smoothness of the coefficients of the un¬ 
derlying equation. Section [4] largely concentrates on constant coefficient higher 
order operators, in particular, the polyharmonic equation, in domains with Lips- 
chitz boundaries. Large parts of this section are taken verbatim from our earlier 
survey [25]; we have added some recent results of Brewster, I. Mitrea and M. Mitrea. 
We have chosen to keep our description of the older results, for completeness, and 
also to provide background and motivation for study of the boundary problems for 
operators with variable coefficients—the subject of Section 0 

The simplest example is the L p -Dirichlet problem for the bilaplacian 

(1.3) A 2 u = 0 in fl, «| an = / € W?(dSl), d„u\ gQ = g e L p (dn), 
in which case the expected sharp estimate on the solution is 

(1.4) l|iV(V W )|| tP(aQ) < C||V r /|| 

LP(d n) + C\\g\\ 

ip(an)’ 

where N denotes the non-tangential maximal function and W[(dfi) is the Sobolev 
space of functions with one tangential derivative in L p (cf. Section [2] for precise 
definitions). In Sections |4. II 14. 6l we discuss (11.311 and (11.411 . and more general higher- 
order homogeneous Dirichlet and regularity boundary value problems with constant 
coefficients, with boundary data in L p . Section 14.71 describes the specific case 
of convex domains. The Neumann problem for the bilaplacian is addressed in 
Section 14.81 In Section 14.91 we discuss inhomogeneous boundary value problems 
(such as problem (11.11) . with ft / 0); for such problems it is natural to consider 
boundary data /, g in Besov spaces, which, in a sense, are intermediate between 
those with Dirichlet and regularity data. 

Finally, in Section [5] we discuss higher order operators with non-smooth coeffi¬ 
cients. The results are still very scarce, but the developments of past several years 
promise to lay a foundation for a general theory. 

To begin, let us mention that contrary to the second order scenario, there are 
several natural generalizations of higher order differential equations to the variable 
coefficient context. Recall that the prototypical higher order operator is the bihar¬ 
monic operator A 2 ; there are two natural ways of writing the biharmonic operator, 
either as a composition A 2 u = A(A u), or in higher-order divergence form 

n n „ 

A 2 u = '^2'^2d j d k (d j d k u) = ^ —d a {d a u). 

1=1 fc=1 M =2 

If we regard A 2 as a composition of two copies of the Laplacian, then one generaliza¬ 
tion to variable coefficients is to replace each copy by a more general second-order 
variable coefficient operator = — div A\7 for some matrix A; this yields operators 
in composition form 

(1.5) Lu{X) = div B(X) V(a(A) div A(X)Vu{X)) 

for some scalar-valued function a and two matrices A and B. Conversely, if we 
regard A 2 as a divergence-form operator, we may generalize to a variable-coefficient 
operator in divergence form 

(1.6) Lu(X) = (-l) m div m A\7 m = (—l) m d a {a aP {X)d p u{X)). 

M = l/3|= m 
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Both classes of operators will be defined more precisely in Section [2] we will see that 
the composition form is closely connected to changes of variables, while operators 
in divergence form are directly associated to positive bilinear forms. 

Sections l5.1l and l5.2l discuss higher order operators in divergence form m : these 
sections discuss, respectively, the Kato problem and the known well-posedness re¬ 
sults for higher-order operators, all of which at present require boundary data in 
fractional smoothness spaces (cf. Section l4~9l) . Section IP1 addresses well-posedness 
of the Diriclrlet boundary value problem for a fourth order operator in a composi¬ 
tion form CE3 with data in L p , p = 2, in particular, generalizing the corresponding 
results for the biharmonic problem (USD- To date, this is the only result address¬ 
ing well-posedness with L p boundary data for higher order elliptic operators with 
bounded measurable coefficients, and its extensions to more general operators, other 
values of p, and other types of boundary data are still open. Returning to diver¬ 
gence form operators & one is bound to start with the very foundations of 
the theory—the estimates on the fundamental solutions. This is a subject of Sec¬ 
tion [531 and the emerging results are new even in the second order case. Having 
those at hand, and relying on the Kato problem solution [16] , we plan to pass to the 
study of the corresponding layer potentials and eventually, to the well-posedness 
of boundary value problems with data in L p . In this context, an interesting new 
challenge, unparalleled in the second order case, is a proper definition of “natural” 
Neumann boundary data. For higher order equations the choice of Neumann data 
is not unique. Depending on peculiarities of the Neumann operator, one can be 
led to well-posed and ill-posed problems even for the bilaplacian, and more general 
operators give rise to new issues related to the coercivity of the underlying form. 
We extensively discuss these issues in the body of the paper and present a certain 
functional analytic approach to definitions in Section [5.51 Finally, Section [5.61 lavs 
out open questions and some preliminary results. 

To finish the introduction, let us point out a few directions of analysis of higher 
order operators on non-smooth domains not covered in this survey. First, an excel¬ 
lent expository paper [50] by Vladimir Maz’ya on the topic of the Wiener criterion 
and pointwise estimates details the state of the art in the end of the previous century 
and provides a considerably more extended discussion of the questions we raised 
in Section [3] surrounding results and open problems. Here, we have concentrated 
on recent developments for the polyharmonic equation and their historical context. 
Secondly, we do not touch upon the methods and results of the part of elliptic 
theory studying the behavior of solutions in the domains with isolated singulari¬ 
ties, conical points, cuspidal points, etc. For a good first exposure to that theory, 
one can consult, e.g., m and references therein. Instead, we have intentionally 
concentrated on the case of Lipschitz domains, which can display accumulating 
singularities—a feature drastically affecting both the available techniques and the 
actual properties of solutions. 

2. Higher order operators: divergence form and composition form 

As we pointed out in the introduction, the prototypical higher-order elliptic 
equation is the biharmonic equation A 2 u = 0, or, more generally, the polyhar¬ 
monic equation A m u = 0 for some integer m > 2. It naturally arises in numerous 
applications in physics and in engineering, and in mathematics it is a basic model for 
a higher-order partial differential equation. For second-order differential equations, 
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the natural generalization of the Laplacian is a divergence-form elliptic operator. 
However, it turns out that even defining a suitable general higher-order elliptic 
operator with variable coefficients is already a challenging problem with multiple 
different solutions, each of them important in its own right. 

Recall that there are two important features possessed by the polyharmonic 
operator. First, it is a “divergence form” operator in the sense that there is an 
associated positive bilinear form, and this positive bilinear form can be used in a 
number of ways; in particular, it allows us to define weak solutions in appropri¬ 
ate Sobolev space. Secondly, it is a “composition operator”, that is, it is defined 
by composition of several copies of the Laplacian. Moreover, if one considers the 
differential equation obtained from the polyharmonic equation by change of vari¬ 
ables, the result would again be a composition of second-order operators. Hence, 
both generalizations are interesting and important for applications, albeit leading 
to different higher-order differential equations. 

Let us discuss the details. To start, a general constant coefficient elliptic operator 
is defined as follows. 

Definition 2.1. Let L be an operator acting on functions u : R" i —> C f . Suppose 
that we may write 

t 

(2.2) (Lu)j = E E 9 a a^ p d 0 u k 

k—1 \a\ = \{3\=m 

for some coefficients a)5 defined for all 1 < j,k < £ and all multiindices a , j3 of 
length n with \a\ = |/3| = m. Then we say that L is a differential operator of 
order 2m. 

Suppose the coefficients are constant and satisfy the Legendre-Hadamard 
ellipticity condition 

(2.3) Re E 53 > A|£| 2m |C| 2 

j,k=l \a\ = \(3\=m 

for all £ £ R" and all £ £ C^, where A > 0 is a real constant. Then we say that L 
is an elliptic operator of order 2m. 

If £ = 1 we say that L is a scalar operator and refer to the equation Lu = 0 as 
an elliptic equation; if £ > 1 we refer to Lu = 0 as an elliptic system. If a jfc = a k f 
then we say the operator L is symmetric. If (f k p is real for all a, (3, j, and k , we 
say that L has real coefficients. 

Here if a is a multiindex of length n, then d a = dff ... 9“^. 

Now let us discuss the case of variable coefficients. A divergence-form higher- 
order elliptic operator is given by 

t 

(2.4) (Lu) j (X) = J 3 53 9>^(A')9VP0). 

k—1 \a\ = \(3\=m 

If the coefficients : R n —> C are sufficiently smooth, we may rewrite m in 
nondivergence form 

t 

(Lu)j(X) = 53 E < k (X)d a u k (X). 

k—1 \a\<2m 


(2.5) 
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This form is particularly convenient when we allow equations with lower-order terms 
(note their appearance in (12.51) 1. 

A simple criterion for ellipticity of the operators L of (12.51) is the condition that 
(12.31) holds with replaced by a^(X) for any X € 1", that is, that 


t 

(2.6) Re E E > Aici 2m ici 2 

j,k— 1 \a\—2m 

for any fixed X £ R ra and for all £ £ R", £ £ C 2 . This means in particular that 
ellipticity is only a property of the highest-order terms of (12.51) : the value of ajf 5 , 
for |a| < m, is not considered. 

Returning to divergence form operators (EH), notice that in this case we have a 
notion of weak solution; we say that Lu = h weakly if 


(2.7) 



E Ec-W 

|a| = |/3|=mj,fc=l JQ 


for any function ip : 0 i— > C f smooth and compactly supported. The right-hand 
side (</?, Lu) may be regarded as a bilinear form A[ip, u] . If L satisfies the ellipticity 
condition (12.6p . then this bilinear form is positive definite. A more general ellipticity 
condition available in the divergence case is simply that (ip,Ltp) > A||V m <£|| i a for 
all appropriate test functions ip (see formula (15.31) below); this condition is precisely 
that the form A[<p, it] be positive definite. 

As mentioned above, there is another important form of higher-order operators. 
Observe that second-order divergence-form equations arise from a change of vari¬ 
ables as follows. If Ait = /i, and u = u o p for some change of variables p , then 

a div AVu = h , 


where a(X) is a real number and A(X) is a real symmetric matrix (both depend¬ 
ing only on p; see Figure |T|). In particular, if u is harmonic then u satisfies the 
divergence-form equation 

E d a {a a p{X)d p u{X)) = 0 
M=l/3|=i 

and so the study of divergence-form equations in simple domains (such as the upper 
half-space) encompasses the study of harmonic functions in more complicated, not 
necessarily smooth, domains (such as the domain above a Lipschitz graph). 

If A 2 it = 0, however, then after a change of variables u does not satisfy a 
divergence-form equation (that is, an equation of the form (12.41) 1. Instead, u satisfies 
an equation of the following composition form: 

(2.8) div AV(adiv AVu) = 0. 

In Section IQ1 we shall discuss some new results pertaining to such operators. 

Finally, let us mention that throughout we let C and e denote positive constants 
whose value may change from line to line. We let f denote the average integral, that 
^ s ’ / dl Jj = -p^E) Ie f onl y measures we will consider are the Lebesgue 

measure dX (on R” or on domains in M") or the surface measure da (on the 
boundaries of domains). 
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Figure 1. The behavior of Laplace’s equation after change of vari¬ 
ables. Here u = uop and J p is the Jacobean matrix for the change 
of variables p. 


3. Boundedness and continuity of derivatives of solutions, the 

MAXIMUM PRINCIPLE AND THE WlENER TEST 


3.1. Miranda-Agmon maximum principle and related geometric restric¬ 
tions on the boundary. The maximum principle for harmonic functions is one 
of the fundamental results in the theory of elliptic equations. It holds in arbi¬ 
trary domains and guarantees that every solution to the Dirichlet problem for the 
Laplace equation, with bounded data, is bounded. Moreover, it remains valid for 
all second-order divergence-form elliptic equations with real coefficients. 

In the case of equations of higher order, the maximum principle has been es¬ 
tablished only in relatively nice domains. It was proven to hold for operators with 
smooth coefficients in smooth domains of dimension two in m and 11041 . and of 
arbitrary dimension in [BJ. In the early 1990s, it was extended to three-dimensional 
domains diffeomorphic to a polyhedron ([73) ( lOOj ) or having a Lipschitz boundary 
f |118l IT201 1. However, in general domains, no direct analog of the maximum prin¬ 
ciple exists (see Problem 4.3, p. 275, in Necas’s book 1140 . The increase of the 
order leads to the failure of the methods which work for second order equations, 
and the properties of the solutions themselves become more involved. 

To be more specific, the following theorem was proved by Agmon. 


Theorem 3.1 ([6} Theorem 1]). Let m > 1 be an integer. Suppose that 12 is domain 
with C 2m boundary. Let 

L= a a {X)d a 

\cx.\<2m 

be a scalar operator of order 2m, where a a £ Cl“l(12). Suppose that L is elliptic 
in the sense of m ■ Suppose further that solutions to the Dirichlet problem for L 
are unique. 

Then, for every u £ C m ~ 1 ( 12) fl C 2rra (12) that satisfies Lu = 0 in 12, we have 
(3.2) max ||<9“u|| ioo(n , < C max ||9^u|| L oo f n nv 

M<m-1 L { ' il > \p\<m—l L 


We remark that the requirement that the Dirichlet problem have unique solu¬ 
tions is not automatically satisfied for elliptic equations with lower-order terms; 
for example, if A is an eigenvalue of the Laplacian then solutions to the Dirichlet 
problem for An — Xu are not unique. 

Equation (13.21) is called the Agmon-Miranda maximum principle. In [1241 , Sul’ce 
generalized this to systems of the form (12.51) , elliptic in the sense of (12.61) , that satisfy 
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a positivity condition (strong enough to imply Agmon’s requirement that solutions 
to the Dirichlet problem be unique). 

Thus the Agmon-Miranda maximum principle holds for sufficiently smooth op¬ 
erators and domains. Moreover, for some operators, the maximum principle is valid 
even in domains with Lipschitz boundary, provided the dimension is small enough. 
We postpone a more detailed discussion of the Lipschitz case to Section l4~6l here we 
simply state the main results. In [1181 and }120| . Pipher and Verchota showed that 
the maximum principle holds for the biharmonic operator A 2 , and more generally 
for the polyharmonic operator A m , in bounded Lipschitz domains in R 2 or R 3 . 
In [141 . Section 8], Verchota extended this to symmetric, strongly elliptic systems 
with real constant coefficients in three-dimensional Lipschitz domains. 

For Laplace’s equation and more general second order elliptic operators, the 
maximum principle continues to hold in arbitrary bounded domains. In contrast, 
the maximum principle for higher-order operators in rough domains generally fails. 

In [96], Maz’ya, Nazarov and Plamenevskii studied the Dirichlet problem (with 
zero boundary data) for constant-coefficient elliptic systems in cones. Counterex¬ 
amples to (13.211 for systems of order 2m in dimension n > 2m+1 immediately follow 
from their results. (See [961 formulas (1.3), (1.18) and (1.28)].) Furthermore, Pipher 
and Verchota constructed counterexamples to (13.21) for the biharmonic operator A 2 
in dimension n = 4 in ( 1171 Section 10], and for the polyharmonic equation A m u = 0 
in dimension n, 4 < n < 2 m+ 1, in [120] Theorem 2.1]. Independently Maz’ya and 
Rossmann showed that (13.21) fails in the exterior of a sufficiently thin cone in di¬ 
mension n, n > 4, where L is any constant-coefficient elliptic scalar operator of 
order 2m > 4 (without lower-order terms). See 11011 Theorem 8 and Remark 3]. 

Moreover, with the exception of [1011 Theorem 8], the aforementioned coun¬ 
terexamples actually provide a stronger negative result than simply the failure of 
the maximum principle: they show that the left-hand side of (13.21) may be infinite 
even if the data of the elliptic problem is as nice as possible, that is, smooth and 
compactly supported. 

The counterexamples, however, pertain to high dimensions. This phenomenon 
raises two fundamental questions: whether the boundedness of the (m— l)-st deriva¬ 
tives remains valid in dimensions n < 3, and whether there are some other, possibly 
lower-order, estimates that characterize the solutions when n > 4. This issue has 
been completely settled in |83l and [85] for the polyharmonic equation in arbitrary 
domains. 

3.2. Sharp pointwise estimates on the derivatives of solutions in arbi¬ 
trary domains. The main results addressing pointwise bounds for solutions to 
the polyharmonic equation in arbitrary domains are as follows. 

Theorem 3.3 ([85]). Let Q be a bounded domain in R", 2 < n < 2 m + 1, and 

(3.4) (-A ) m u = finn, f e C7g°(fi), u £ IV m ’ 2 (fI). 

Then the solution to the boundary value problem HMD satisfies 

(3.5) V m_n / 2+1 / 2 u g L°°(fl) when n is odd, V m-ra//2 u e L°°( fl) when n is even. 
In particular, 

(3.6) V m-1 u e L°°{n) when n = 2,3. 
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Here the space W m,2 (I2), is, as usual, a completion of Cg°(fl) in the norm given 
by = ||V m zi||i 2 (Q). We note that W m,2 (f2) embeds into C k (fl) only 

when k is strictly smaller than m — n < 2m. Thus, whether the dimension 
is even or odd, Theorem 13.31 gains one derivative over the outcome of Sobolev 
embedding. 

The results of Theorem 13.31 are sharp, in the sense that the solutions do not ex¬ 
hibit higher smoothness than warranted by (13.51) (13.61) in general domains. Indeed, 
assume that n £ [3,2 m + 1] D N is odd and let H C R n be the punctured unit ball 
B\ \ {O}, where B r = {x £ R" : |x| < r}. Consider a function 77 £ Cq°(B 1 / 2 ) such 
that 77 = 1 on B\u. Then let 

(3.7) u(x) :=77(*)ar" t_ M|x| 2m - n ), x e B ± \ {O}, 

where d x stands for a derivative in the direction of Xi for some i = 1,... ,n. It 
is straightforward to check that u £ W m ’ 2 (fl) and (—A) m w £ Co°(f2). While 
V m_ v +2 u is bounded, the derivatives of order m — ^ | are not, and moreover, 

V m-: i+5u is not continuous at the origin. Therefore, the estimates (13.51) (13.61) are 
optimal in general domains. 

As for the case when n is even, the results in E3 Section 10.4] demonstrate that 
in the exterior of a ray there is an m-harmonic function behaving as |a;| m_ 2 '+ 2 . 
Thus, upon truncation by the aforementioned cut-off 77, one obtains a solution to 
(13.41) in = 0,..., x n -\ = 0,0 < x n < 1}, whose derivatives of order m—^ + 1 

are not bounded. More delicate examples can be obtained from our results for the 
Wiener test to be discussed in Section ET41 Those show that the derivatives of order 
777. — ^ need not be continuous in even dimensions. Therefore, in even dimensions 
(13.51) is a sharp property as well. 

It is worth noting that the results above address also boundedness of solutions 
(rather than their derivatives) corresponding to the case the case when 777.— j+ i = 0 
in odd dimensions, or, respectively, m — j = 0 in the even case. In this respect, 
we would also like to mention higher dimensional results following from the Green 
function estimates in jJQ]. As will be discussed in the next section, one can show 
that, in addition to our results above, if S2 C R” is bounded for n < 2m + 2, and 
if u is a solution to the polyharmonic equation (13.41) . then u £ L°°(fl). This result 
also holds if 12 C R 7 and m = 2. 

If 12 C R™ is bounded and n > 2m + 3, or if m = 2 and n> 8, then the question 
of whether solutions u to (ED) are bounded is open. In particular, it is not known 
whether solutions u to 

A 2 7t, = h in 12, u £ W’|(f2) 

are bounded if 12 C R" for n > 8. However, there exists another fourth-order 
operator whose solutions are not bounded in higher-dimensional domains. In 541. 
Maz’ya and Nazarov showed that if n > 8 and if a > 0 is large enough, then there 
exists an open cone K C R" and a function h £ Co°(iv \{0}) such that the solution 
u to 

(3.8) A 2 it + ad^u = h in K, u £ W 2 (K) 
is unbounded near the origin. 

3.3. Green function estimates. Theorem 13.31 has several quantitative manifes¬ 
tations, providing specific estimates on the solutions to (13.41) . Most importantly, 
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the authors established sharp pointwise estimates on Green’s function of the poly¬ 
harmonic operator and its derivatives, once again without any restrictions on the 
geometry of the domain. 

To start, let us recall the definition of the fundamental solution for the polyhar¬ 
monic equation (see. e.g., [10] 1. A fundamental solution for the m-Laplacian is a 
linear combination of the characteristic singular solution (defined below) and any 
m- harmonic function in R n . The characteristic singular solution is 

(3.9) C m ,n\x\ 2m ~ n , if n is odd, or if n is even with n > 2m + 2, 

(3.10) C m ,n\x\ 2rn ~ n log |x|, if n is even with n < 2m. 

The exact expressions for constants C mi „ can be found in m, p. 8. Hereafter we 
will use the fundamental solution given by 

{ \x\ 2m ~ n , if n is odd, 

| x\ 2m ~ n log if n is even and n < 2m, 

\x\ 2m -n^ if n is even and n > 2m + 2. 

As is customary, we denote the Green’s function for the polyharmonic equation 
by G(x,y ), x, y € 12, and its regular part by S(x,y), that is, S(x,y ) = G(x,y) — 
T(x — y). By definition, for every fixed y € ft the function G(- ,y) satisfies 

(3.12) (—A x ) m G(x, y) = S(x -y), x £ ft, 

in the space IT™’ 2 (ft). Here A^ stands for the Laplacian in the x variable. Similarly, 
we use the notation A y , V, y , V x for the Laplacian and gradient in y , and gradient 
in x, respectively. By d(x) we denote the distance from a; £ ft to <9ft. 


Theorem 3.13. Let ft C KG be an arbitrary bounded domain, m £ N, n £ [2, 2m + 
1] fl N, and let 

{ m — n/2 + 1/2 when n is odd , 

m — n/2 when n is even. 


Fix any number N > 25. Then there exists a constant C depending only on m, n, 
N such that for every x,y £ ft the following estimates hold. 

If n £ [3, 2m + 1] fl N is odd then 
(3.15) 

d(v')^ ^ 

I KV 3 y G(x, y)\ < C |a; _ y ^ +n _ 2m+i , when\x-y\>Nd(y), 0<i,j <\, 

and 

(3 ' 16) 

\S7 l x S7 J y G(x,y)\ < C _^| n _ 2m+j , when \x - y\ > Nd{x), 0<i,j<\. 

Next, 

P-U) |V‘V»G(X,;,)|< k i<| „ C 2m+ . +J . , 


when | a; 
and 


(3.18) 


y \ < N 1 max{d(a;), d(y)}, and i + j > 2m — n, 0 < i, j < m — n/2 + 1/2, 
V'XfGU.y) < C min{<i(a:), d(y)} 2m ~ n ~ i ~\ 
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when \x — y\ < N 1 max{d(x), d(y)}, and i + j< 2 m — n, 0 < i, j < to — n/2 + 1/2. 
Finally, 


(3.19) 


\V l x V y G(x,y)\< 


C 


mm{d(x), d(y), \x — y\}n-2m+i+j 

__ C _ 

ma x{d(x),d(y), |x — y\}ri-2m+i+j ’ 

when N~ 1 d(x) < \x — y\ < Nd(x) and N~ 1 d{y) < \x — y\ < Nd(y), 0 < i, j < A. 

Furthermore, if n £ [3, 2m + 1] (~l N is odd, the estimates on the regular part of 
the Green function S are as follows: 

(3.20) 

C 

when \x — y\ > N min{d(x),d(y)}, 0 < i,j < A. 


\^ v S{x-y)\ < 
Next, 

(3.21) 


|j; _ y\n—2m+i+j 

|ViVj5(s,i/)| < 


C 


max{d(x), d(y)} n 2m + l +o 


when \x — y\ < N 1 max{<i(:r), d(y)}, and i + j> 2 m — n, 0 < i, j < m — n/2 + 1/2, 
and 


(3.22) 


|V^V^5(x,y)| < Cmm{d(x),d{y)} 


2m—n—i—j 


when \x — y\ < N 1 max{d(a;), d(y)}, and i + j< 2m — n, 0 < i, j < m — n/2 + 1/2. 
Finally, 


(3.23) 


|KV^0r,2/)|< — 


C 


min{d(ir), d(y), \x — y|} n - 2m +*+l 

__ C _ 

max{d(x),d(y), \x — y\}n-2m+i+j ’ 

when N~ x d(x) < \x — y\ < Nd(x) and N^ 1 d(y) < \x — y\ < Nd(y), 0 < i,j < 
m — n/2 + 1/2. 

Ifn £ [2, 2m]flN is even, then (13.1511 (13.1611 and (13.191) are valid with A = to— j, 
and 

(3.24) 

|ViVjG(aM,)| < Cunn{d(x),d(y)} 2m - n -^ (c 1 + log , 

when \x — y\ < max{<i(x), d(y)} and 0 < i, j < m — n/2. 

Furthermore, if n £ [2, 2to] fl N is even, the estimates on the regular part of the 
Green function S are as follows: 

(3.25) |V' x VZS(x -y)l<C\x- y\-^-i-o ^ + log 


when \x — y\ > N min{d(x), d(y)}, 0 < i, j < to — n/2. Next, 
(3.26) 

|Vl.V{S(x, 2 /)| < Cmm{d(x),d(y)} 2m - n - i ^ (c' + log 


diam fl 


max{d(x), d(y)} J 
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when \x — y\ < N 1 max{d(x), d(y)}, 0 < i, j < m — n/2. Finally, 


(3.27) \V' x VlS(x, y)\ < C min {d(x),d(y), \x - y |} 


2 m—n—i—j > 


x [C + log 


diam Q 


ma x{d(x),d(y), \x — y\}n-2m+i+j 

when N~ x d(x) < \x — y\ < Nd(x) and TV -1 d(y) < \x — y\ < Nd(y), 0 < i,j < 
m — n/2. 

We would like to highlight the most important case of the estimates above, 
pertaining to the highest order derivatives. 

Corollary 3.28. Let 12 C R n be an arbitrary bounded domain. If n G [3, 2m+l]nN 
is odd, then for all x, y G LI, 


(3.29) 
and, in particular, 

(3.30) 


V™ * + *(G(x,y)-T(x-y)) 


< 


C 


Vr S+5 G(a;,y) 


< 


max{d(x),d(y),\x - y\} 7 
C 


~y r 


If n G [2, 2m] D N is even, then for all x, y G LI, 


(3.31) 

and 

(3.32) 


-,m— ’rn— 


Vy *[G[x,y)-T{x-y)) 

< C log ( 1 + 


diam H 


max{d(x), d(y), \x — y\} J " 


y m-~ m- 
\ x V y 


^(x, ri |< C iog(i i -wy 1 ). 


The constant C in (13.29P (13.321) depends on m and n only. In particular, it does 
not depend on the size or the geometry of the domain f2. 


We mention that the pointwise bounds on the absolute value of Green’s function 
itself have been treated previously in dimensions 2 m + 1 and 2m + 2 for m > 2 
and dimensions 5, 6, 7 for m = 2 in [901 Section 10] (see also [87!). I n particular, in 
H3 Section 10], Maz’ya showed that the Green’s function G m (x,y) for A m in an 
arbitrary bounded domain f2 C R” satisfies 

(3-33) \G m (x,y)\ < - ^ J_ 2m 

I® ^ V I 

if n = 2m + 1 or n = 2m + 2. If m = 2, then (13.331) also holds in dimension 
n = 7 = 2m + 3 (cf. [87]). Whether (13.331) holds in dimension n > 8 (for m = 2) or 
n > 2m + 3 (for m > 2) is an open problem; see [10 Problem 2], Also, similarly to 
the case of general solutions discussed above, there exist results for Green functions 
in smooth domains @5], 0; [135] . [136] . in conical domains [58], [73], and in 
polyhedra [lOO] , 

Furthermore, using standard techniques, the Green’s function estimates can be 
employed to establish the bounds on the solution to (13.41) for general classes of data 
/, such as L p for a certain range of p , Lorentz spaces etc. A sample statement to 
this effect is as follows. 
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Proposition 3.34. Let Ll C M™ be an arbitrary bounded domain, m G N, n G 
[2,2to+ 1] n N, and let A retain the significance of (|3.14|) . Consider the boundary 
value problem 

(3.35) (-A ) m u = £ c a d a f a , u G W m ’ 2 {Ll). 

I«|<A 


Then the solution satisfies the following estimates. 
If n G [3, 2to +l]flN is odd, then for all x G Ll, 


(3.36) |V"*-t+*u(x)| < C m ,. 


E l 

' n . l " is 


d(y) 


m — -f+l — |a| 


\x - 2/| 


I fa (//) I dy, 


\a\<m— ■j + l' ' 

whenever the integrals on the right-hand side of (13.361) are finite. In particular, 


(3.37) ||V"-5+* U || £ « (n) 


<C m , n ,n £ ||d(.r-®-i- |a| /a||LP( n)) P> 

|a|<m-f+ i 

provided that the norms on the right-hand side of (13.371) are finite. 

If n G [2, 2to] D N is even, i/ien for all x G LI, 

(3.38) | it ( 2 ;) | 


< C m , n £ J H log (l + |/a(2/)| dy, 

whenever the integrals on the right-hand side of (|3.38[) are finite. In particular, 


(3.39) ||V m -t u || LOO(a) <C m , n , n £ ||d(.)"*-*-l a| / a || L , ( n), p > 1, 

provided that the norms on the right-hand side of (13.391) are finite. 

The constants C m(n above depend on m and n only, while the constants denoted 
by Cm^n^n depend on m, n, and the diameter of the domain LI. 


By the same token, if (13.331) holds, then solutions to (13.41) satisfy 

IMI L oo (n) < C (m, n,p) diam(f]) 2m_ "/ p ||/1| LP(an) 

provided p > n/2m (see, e.g., j5D] Section 2]). Thus, e.g., if Ll C R" is bounded 
for n = 2m + 2, and if u satisfies m for a reasonably nice function /, then 
u G L°°(Lt). This result also holds if Ll C K 7 and m = 2. This complements the 
results in Theorem 13.31 as discussed in Section m 

To conclude our discussion of Green’s functions, we mention two results from 
[105] ; these results are restricted to relatively well-behaved domains. In [TO 5] . 
D. Mitrea and I. Mitrea showed that, if Ll is a bounded Lipschitz domain in R 3 , 
and G denotes the Green’s function for the bilaplacian A 2 , then the estimates 

V 2 G(z, •) e L 3 (Ll), dist( •, dLl)~ a WG{x, ■) e L 3 /“’°° 

hold, uniformly in x G Ll, for all 0 < a < 1. 

Moreover, they considered more general elliptic systems. Suppose that L is 
an arbitrary elliptic operator of order 2m with constant coefficients, as defined 
by Definition 12.11 and that G denotes the Green’s function for L. Suppose that 
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Cl C M™, for n > m, is a Lipscliitz domain, and that the unit outward normal v to 
f l lies in the Sarason space VA'IO(dCl) of functions of vanishing mean oscillations 
on dCl. Then the estimates 

(3.40) V m G(x, ■) G -°°(n), 

dist(•,5fi)- a V m - 1 G'( x, •) G £ - m -i+. .°°(n) 
hold, uniformly in x G Cl, for any 0 < a < 1. 

3.4. The Wiener test: continuity of solutions. In this section, we discuss 
conditions that ensure that solutions (or appropriate gradients of solutions) are 
continuous up to the boundary. These conditions parallel the famous result of 
Wiener, who in 1924 formulated a criterion that ensured continuity of harmonic 
functions at boundary points |144j . Wiener’s criterion has been extended to a 
variety of second-order elliptic and parabolic equations ((77] [50] 09j 021 [79l [3] 11381 
(7511481 : see also the review papers mm- However, as with the maximum principle, 
extending this criterion to higher-order elliptic equations is a subtle matter, and 
many open questions remain. 

We begin by stating the classical Wiener criterion for the Laplacian. If C R" 
is a domain and Q G dCl, then Q is called regular for the Laplacian if every solution 
u to 

A u = h in Cl, u G W 2 (Cl) 

for h G Cq°(CI) satisfies lim x^-qu(X) = 0. According to Wiener’s theorem [144j . 
the boundary point Q G dCl is regular if and only if the equation 

(3.41) f cap 2 (.B(Q, s) \ fl)s 1_ri ds = oo 
Jo 

holds, where 

(3.42) cap 2 (/v) = inf{|M|£ 2(Hn) + ||Vu||^ 2(Rn) : u G C^°(R n ), u > 1 on if}. 

For example, suppose Cl satisfies the exterior cone condition at Q. That is, 
suppose there is some open cone K with vertex at Q and some e > 0 such that 
if fl B(Q 7 s) C Cl c . It is elementary to show that cap 2 (B(Q, s) \Cl) > C(K)s n ~ 2 
for all 0 < s < £, and so (13.4111 holds and Q is regular. Regularity of such points 
was known prior to Wiener (see [121 ,. [1451 . and [76]) and provided inspiration for 
the formulation of the Wiener test. 

By m if L = — div AV is a second-order divergence-form operator, where the 
matrix A(X) is bounded, measurable, real, symmetric and elliptic, then Q G dCl 
is regular for L if and only if Q and Cl satisfy (13.4111 . In other words, Q G dCl is 
regular for the Laplacian if and only if it is regular for all such operators. Similar 
results hold for some other classes of second-order equations; see, for example, m, 

[32, or [451 . 

One would like to consider the Wiener criterion for higher-order elliptic equa¬ 
tions, and that immediately gives rise to the question of natural generalization of 
the concept of a regular point. The Wiener criterion for the second order PDEs 
ensures, in particular, that weak W 2 solutions are classical. That is, the solution 
approaches its boundary values in the pointwise sense (continuously). From that 
point of view, one would extend the concept of regularity of a boundary point as 
continuity of derivatives of order to — 1 of the solution to an equation of order 2 to up 
to the boundary. On the other hand, as we discussed in the previous section, even 
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the boundedness of solutions cannot be guaranteed in general, and thus, in some 
dimensions the study of the continuity up to the boundary for solutions themselves 
is also very natural. We begin with the latter question, as it is better understood. 

Let us first define a regular point for an arbitrary differential operator L of order 
2m analogously to the case of the Laplacian, by requiring that every solution u to 

(3.43) Lu = h in f2, u £ 

for h £ C^°(f2) satisfy limx~>Qu(X) = 0. Note that by the Sobolev embedding 
theorem, if C R n for n < 2m — 1, then every u £ W^( Cl) is Holder continuous 
on and so satisfies liniA'->.Q u(X) = 0 at every point Q £ dfl. Thus, we are only 
interested in continuity of the solutions at the boundary when n > 2m. 

In this context, the appropriate concept of capacity is the potential-theoretic 
Riesz capacity of order 2m, given by 

(3.44) cap 2m (A') = inf j ^ II^II^Rn) : u e C 0 “(r), u > 1 on Ar|. 

0< |cn| <m 

The following is known. If m > 3, and if C R" for n = 2m, 2m + 1 or 2m + 2, 
or if m = 2 and n = 4, 5, 6 or 7, then Q £ dfl is regular for A m if and only if 

(3.45) f cap 2rn (I?((3, s) \ ds = oo. 

Jo 

The biharmonic case was treated in [86] and [87], and the polyharmonic case for 
m > 3 in [92] and [89]. 

Let us briefly discuss the method of the proof in order to explain the restrictions 
on the dimension. Let L be an arbitrary elliptic operator, and let F be the funda¬ 
mental solution for L in R” with pole at Q. We say that L is positive with weight 
F if, for all u £ C“(R n \ {Q}), we have that 

P m p 

(3.46) / Lu(X)-u(X)F(X)dX>cY J \X k u(X)\ 2 \X\ 2k ~ n dX. 

J\ R n j -j J M. n 

The biharmonic operator is positive with weight F in dimension n if 4 < n < 7, 
and the polyharmonic operator A m , m > 3, is positive with weight F in dimension 
2m < n < 2m -I- 2. (The Laplacian A is positive with weight F in any dimension.) 
The biharmonic operator A 2 is not positive with weight F in dimensions n > 8, and 
A m is not positive with weight F in dimension n > 2m + 3. See [HEi Propositions 
1 and 2]. 

The proof of the Wiener criterion for the polyharmonic operator required pos¬ 
itivity with weight F. In fact, it turns out that positivity with weight F suffices 
to provide a Wiener criterion for an arbitrary scalar elliptic operator with constant 
coefficients. 

Theorem 3.47 ( [91], Theorems 1 and 2]). Suppose Q C R™ and that L is a scalar 
elliptic operator of order 2m with constant real coefficients, as defined by Defini- 
tion \2.1\ 

If n = 2m, then Q £ d£l is regular for L if and only if (13.45[) holds. 

If n > 2m + 1, and if the condition (13.461) holds, then again Q £ dfl is regular 
for L if and only if (13.451) holds. 

This theorem is also valid for certain variable-coefficient operators in divergence 
form; see the remark at the end of [89] Section 5]. 
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Similar results have been proven for some second-order elliptic systems. In par¬ 
ticular, for the Lame system Lu = Au+agraddivu, a > — 1, positivity with weight 
F and Wiener criterion have been established for a range of a close to zero, that 
is, when the underlying operator is close to the Laplacian ([78]). It was also shown 
that positivity with weight F may in general fail for the Lame system. Since the 
present review is restricted to the higher order operators, we shall not elaborate on 
this point and instead refer the reader to m for more detailed discussion. 

In the absense of the positivity condition (13.461) . the situation is much more 
involved. Let us point out first that the condition (13.461) is not necessary for reg¬ 
ularity of a boundary point, that is, the continuity of the solutions. There exist 
fourth-order elliptic operators that are not positive with weight F whose solutions 
exhibit nice behavior near the boundary; there exist other such operators whose 
solutions exhibit very bad behavior near the boundary. 

Specifically, recall that (13.461) fails for L = A 2 in dimension n > 8. Nonetheless, 
solutions to A 2 u = h are often well-behaved near the boundary. By [99] . the 
vertex of a cone is regular for the bilaplacian in any dimension. Furthermore, if the 
capacity condition (13.451) holds with m = 2, then by [91] Section 10], any solution 
u to 

A 2 u = h in 9, u £ (9) 

for h £ C'o°(f2) satisfies limx^Q u(X) = 0 provided the limit is taken along a 
nontangential direction. 

Conversely, if n > 8 and L = A 2 +ad then by [94] . there exists a cone K and a 
function h £ Cq°(K\{0}) such that the solution u to (13.81) is not only discontinuous 
but unbounded near the vertex of the cone. We remark that a careful examination 
of the proof in [53] implies that solutions to (13.81) are unbounded even along some 
nontangential directions. 

Thus, conical points in dimension eight are regular for the bilaplacian and irregu¬ 
lar for the operator A 2 +a c^. Hence, a relevant Wiener condition must use different 
capacities for these two operators. This is a striking contrast with the second-order 
case, where the same capacity condition implies regularity for all divergence-form 
operators, even with variable coefficients. 

This concludes the discussion of regularity in terms of continuity of the solution. 
We now turn to regularity in terms of continuity of the (m — l)-st derivatives. 
Unfortunately, much less is known in this case. 


3.5. The higher order Wiener test: continuity of derivatives of poly- 
harmonic functions. The most natural generalization of the Wiener test to the 
higher order scenario concerns the continuity of the derivatives of the solutions, 
rather than solutions themselves, as derivatives constitute part of the boundary 
data. However, a necessary prerequisite for such results is boundedness of the cor¬ 
responding derivatives of the solutions—an extremely delicate matter in its own 
right as detailed in Section 13.21 In the context of the polyharmonic equation, The- 
orem limi has set the stage for an extensive investigation of the Wiener criterion and, 
following earlier results in [82], the second author of this paper and V. Maz’ya have 
recently obtained a full extension of the Wiener test to the polyharmonic context in 
[531 . One of the most intricate issues is the proper definition of the polyharmonic 
capacity, and we start by addressing it. 
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At this point Theorem 13.31 finally sets the stage for a discussion of the Wiener 
test for continuity of the corresponding derivatives of the solution, which brings us 
to the main results of the present paper. 

Assume that m?N and n £ [2, 2m + 1] D N. Let us denote by Z the following 
set of indices: 

(3.48) Z ={0,1,... ,m — n/2 + 1/2} if n is odd, 

(3.49) Z ={-n/ 2 + 2, - n/2 + 4,..., m - n/2 - 2, m - n/2} n (N U {0}) 

if n is even, m is even, 

(3.50) Z ={—n/2 + 1, —n/2 + 3,..., m — n/2 — 2, to — n/2}fl(NU {0}) 

if n is even, m is odd. 

Now let II be the space of linear combinations of spherical harmonics 

v 

(3.51) P(x) = EEv Y*(x/ |x|), b pl eR, xeR n \ {O}, 

pGZ l=-p 

with the norm 

(3.52) nin:=(EE^) and 11, := {P € II: ||P|| n = 1}. 

\p<EZl=-p J 

Then, given P £ III, an open set D in R™ such that O £ R™\D, and a compactum 
K in D, we define 

(3.53) Cap P (K, D) 


:= inf | J |V m n(x)| 2 dx : u £ W m,2 (D), u = P in a neighborhood of K 
with 

(3.54) Cap (K,D):= inf Cap P (K,D). 

PGlli 

In the context of the Wiener test, we will be working extensively with the capacity 
of the complement of a domain Q C R” in the balls B 2 -j , j £ N, and even more 
so, in dyadic annuli, C 2 -i i 2 -i+ 2 , j £ N, where C SiQS := {x £ R" : s < |x| < as}, s, 
a > 0. As is customary, we will drop the reference to the “ambient” set 

(3.55) Capp (C 2 -i, 2 - 1+2 \ Cl) := Capp (C 2 -j t2 -i+ 2 \H, C 2 -i-2 2 - 3 +4), j £ N, 

and will drop the similar reference for Cap. In fact, it will be proven below that there 
are several equivalent definitions of capacity, in particular, for any n £ [2, 2m + 1] 
and for any s > 0, a > 0, K C C Stas , we have 

(3.56) Cap P (K,C s/2i2as ) 

» i" f { E / ESr *■ : « 6 \ (OR, 

U=o jRn |x| 

u = P in a neighborhood of K 
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In the case when the dimension is odd, also 

Capp (C Sj a S \ ft, C s /2,2as) ~ Capp ( C St a S \ ft, R" \ {O}). 

Thus, either of the above can be used in (13.551) . as convenient. 

Let ft be a domain in R ra , n > 2. The point Q £ 9ft is k-regular with respect 
to the domain ft and the operator (—A) m , m £ N, if the solution to the boundary 
problem 

(3.57) (—A) m u = / in ft, / e C 0 °°(ft), u £ W m ’ 2 (ft), 
satisfies the condition 

(3.58) S7 k u(x) — > 0 as x —» Q, x £ ft, 

that is, all partial derivatives of u of order k are continuous. Otherwise, we say 
that Q £ 9ft is /e-irregular. 


Theorem 3.59 ([82], [84]). Let ft be an arbitrary open set in R”, m £ N, 2 < n < 
2to + 1. Let X be given by 


(3.60) 

If 


X = 


m — n/2 + 1/2 when n is odd, 


m — n/2 when n is even. 


OO 


(3.61) 



3=0 

and 


(3.62) 

OO 

j 2-j( 2m - n ) 


3=0 


inf 

Peni 


Capp ((^ 2 —j, 2 — j+ 2 \ ^) — +oo, 


inf 

Peni 


Capp (02-3^2-1+2 


\ ft) = +oo, 


when n is odd, 


when n is even, 


then the point O is X-regular with respect to the domain ft and the operator (—A) m . 

Conversely, if the point O £ 9ft is X-regular with inspect to the domain ft and 
the operator (—A) m then 


(3.63) inf 2 2m Cap P ( 02 - 3 , 2 - 1+2 \ ft) = + 00 , when n is odd, 


(3.64) inf 2 ^ 2m Cap P ( 02 - 1 , 2 - 1+2 \ ft) = +c«, when n is even. 

£ ^'=0 

Here, as before, 0 2 -i 2 - 1+2 is the annulus {x £ R" : 2 3 < |x| < 2 ■ J+2 } ) j £ 
NU{0}. 


Let us now discuss the results of Tlieorem l3.59l in more detail. This was the first 
treatment of the continuity of derivatives of an elliptic equation of order m > 2 
at the boundary, and the first time the capacity (13.531) appeared in the literature. 
When applied to the case m = 1, n = 3, it yields the classical Wiener criterion 
for continuity of a harmonic function (cf. (13.411) 1. Furthermore, as discussed in 
the previous section, continuity of the solution itself (rather than its derivatives) 
has been previously treated for the polyharmonic equation, and for (—A)” 1 the 
resulting criterion also follows from Theorem 13.591 in particular, when m = 2n, the 
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new notion of capacity (13.4811 (13.5211 coincides with the potential-theoretical Bessel 
capacity used in m- In the case A = 0, covering both of the above, necessary 
and sufficient condition in Theorem 13.591 are trivially the same, as P = 1 when 
n = 2m in even dimensions and n = 2m + 1 in odd ones. For lower dimensions 
n the discrepancy is not artificial, for, e.g., (13.6111 may fail to be necessary as was 
shown in [82] . 

It is not difficult to verify that we also recover known bounds in Lipschitz and 
in smooth domains, as the capacity of a cone and hence, capacity of an intersec¬ 
tion with a complement of a Lipschitz domains, assures divergence of the series 
in (13.611) (13.621) . On the other hand, given Theorem 13.591 and following considera¬ 
tions traditional in this context (choosing sufficiently small balls in the consecutive 
annuli to constitute a complement of the domain), we can build a set with a conver¬ 
gent capacitory integral and, respectively, an irregular solution with discontinuous 
derivatives of order A at the point O. Note that this yields further sharpness of 
the results of Theorem l3.31 In particular, in even dimensions, it is a stronger coun¬ 
terexample than that of a continuum (not only m — n/2 + 1 derivatives are not 
bounded, but m — n/2 derivatives might be discontinuous). We refer the reader 
back to Section [3~2l for more details. 

One of the most difficult aspects of proof of Theorem 13.591 is finding a correct 
notion of polyharmonic capacity and understanding its key properties. A peculiar 
choice of linear combinations of spherical harmonics (see (13.481) (13.5011 and (13.511) ') 
is crucial at several stages of the argument, specific to the problem at hand, and 
no alterations would lead to reasonable necessary and sufficient conditions. At the 
same time, the new capacity and the notion of higher-order regularity sometimes 
exhibit surprising properties, such as for instance sensitivity to the affine changes of 
coordinates [82], or the aforementioned fact that in sharp contrast with the second 
order case m one does not expect the same geometric conditions to be responsible 
for regularity of solutions to all higher order elliptic equations. 

It is interesting to point out that despite fairly involved definitions, capacitory 
conditions may reduce to a simple and concise criterion, e.g., in a case of a graph. 
To be precise, let C R 3 be a domain whose boundary is the graph of a function ip, 
and let uj be its modulus of continuity. If 


(3.65) 


f 1 tdt 


then every solution to the biharmonic equation satisfies Vu £ C(H). Conversely, 
for every w such that the integral in (13.651) is convergent, there exists a C 0,u> domain 
and a solution u of the biharmonic equation such that Vu ^ C(Q). In particular, as 
expected, the gradient of a solution to the biharmonic equation is always bounded 
in Lipschitz domains and is not necessarily bounded in a Holder domain. Moreover, 
one can deduce from (13.651) that the gradient of a solution is always bounded, e.g., 
in a domain with w(t) ~ tlog 1,/2 f, which is not Lipschitz, and might fail to be 
bounded in a domain with ui(t ) ss tlogt. More properties of the new capacity and 
examples can be found in [82j, [84] . 





























HIGHER-ORDER ELLIPTIC EQUATIONS 


21 


4. Boundary value problems in Lipschitz domains for elliptic 

OPERATORS WITH CONSTANT COEFFICIENTS 

The maximum principle (13.21) provides estimates on solutions whose boundary 
data lies in L°° . Recall that for second-order partial differential equations with real 
coefficients, the maximum principle is valid in arbitrary bounded domains. The 
corresponding sharp estimates for boundary data in L P 1 1 < p < oo, are much more 
delicate. They are not valid in arbitrary domains, even for harmonic functions, 
and they depend in a delicate way on the geometry of the boundary. At present, 
boundary-value problems for the Laplacian and for general real symmetric elliptic 
operators of the second order are fairly well understood on Lipschitz domains. See, 
in particular, [551 . 

We consider biharmonic functions and more general higher-order elliptic equa¬ 
tions. The question of estimates on biharmonic functions with data in L p was raised 
by Riviere in the 1970s ([2H]), and later Kenig redirected it towards Lipschitz do¬ 
mains in [ 551 l66j . The sharp range of well-posedness in L p , even for biharmonic 
functions, remains an open problem (see [66j Problem 3.2.30]). In this section we 
shall review the current state of the art in the subject, the main techniques that 
have been successfully implemented, and their limitations in the higher-order case. 

Most of the results we will discuss are valid in Lipschitz domains, defined as 
follows. 


Definition 4.1. A domain Q C R” is called a Lipschitz domain if, for every 
Q £ dCl , there is a number r > 0, a Lipschitz function ip : R"^ 1 i—► R with 
||V¥>|| Loo < Af, and a rectangular coordinate system for R" such that 

B(Q , r)fl!l = {(x, s) : x £ l"' 1 , s £ R, |(x, s) — Q)\ < r, and s > p(x)}. 

If we may take the functions p to be C k (that is, to possess k continuous deriva¬ 
tives), we say that ft is a C k domain. 

The outward normal vector to O will be denoted v. The surface measure will be 
denoted a, and the tangential derivative along dQ will be denoted V T . 


In this paper, we will assume that all domains under consideration have con¬ 
nected boundary. Furthermore, if dLl is unbounded, we assume that there is a sin¬ 
gle Lipschitz function tp and coordinate system that satisfies the conditions given 
above; that is, we assume that Q is the domain above (in some coordinate system) 
the graph of a Lipschitz function. 

In order to properly state boundary-value problems on Lipschitz domains, we 
will need the notions of non-tangential convergence and non-tangential maximal 
function. 

In this and subsequent sections we say that u| aQ = / if / is the nontangential 
limit of u, that is, if 


lim 

xaQ, xer(Q) 


u(X) = f(Q) 


for almost every (da) Q £ d 12, where T(Q) is the nontangential cone 


(4.2) T(Q) = {Y £ ft : dist(y, dCl ) < (1 + a)\X - Yj}. 


Here a > 0 is a positive parameter; the exact value of a is usually irrelevant to 
applications. The nontangential maximal function is given by 

(4.3) NF(Q) = sup{|F 1 (A')| : X £ T(Q)}. 
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The normal derivative of u of order to is defined as 

w«= E HQr^uiQ), 

i E, a! 

where d a u(Q) is taken in the sense of nontangential limits as usual. 


4.1. The Dirichlet problem: definitions, layer potentials, and some well- 
posedness results. We say that the L p -Diric.hlet problem for the biharmonic op¬ 
erator A 2 in a domain f1 is well-posed if there exists a constant C > 0 such that, 
for every / £ Wf(dfl) and every g £ L p (dfl), there exists a unique function u that 
satisfies 


(4.4) 


A 2 u = 0 
u = f 
d v u = g 

, l|JV(V«)|| i p (a0 ) < C|| ff || LP{dn) 


in f2, 
on 9f2, 
on dfl, 


The L p -Dirichlet problem for the polyharmonic operator A m is somewhat more 
involved, because the notion of boundary data is necessarily more subtle. We say 
that the Z/'-Dirichlet problem for A m in a domain fl is well-posed if there exists a 
constant C > 0 such that, for every g £ L p (dfl) and every / in the Whitney-Sobolev 
space WA p n _ 1 (dQ). there exists a unique function u that satisfies 


(4.5) 


< 


A m u = 0 

^““lan = f< 

d^u = g 


in f2, 

for all 0 < \a\ < m — 2, 
on Oil, 


\LP(an) 


<C\\g\\ 

LP(dQ. + C E ll^' r /“llip(an)- 

\a\=m —2 


The space WA p n (diV) is defined as follows. 


Definition 4.6. Suppose that f 1 C R" is a Lipschitz domain, and consider arrays 
of functions f = {f a ■ |or| < m — 1} indexed by multiindices a of length n, where 
f a : dfl i-A C. We let WA p n ( dQ ) be the completion of the set of arrays ip = {d a ip : 

|a| < to — 1}, for ip £ Cq°(R"), under the norm 

(4-7) + E ll^' r< ^ a V'lliP(dn)- 

|a|<m—1 |a|—m—1 

If we prescribe d a u = f a on dtt for some / £ WA^dCl), then we are prescribing 
the values of u, Vu,..., V m-1 w on dfl, and requiring that (the prescribed part of) 
V m u\ gQ lie in L p (dfl). 

The study of these problems began with biharmonic functions in C 1 domains. 
In [126] , Selvaggi and Sisto proved that, if fl is the domain above the graph of a 
compactly supported C 1 function ip, with ||V</>|| ioo small enough, then solutions 
to the Dirichlet problem exist provided 1 < p < oo. Their method used certain 
biharmonic layer potentials composed with the Riesz transforms. 

In [34], Cohen and Gosselin proved that, if fl is a bounded, simply connected C 1 
domain contained in the plane R 2 , then the L p -Dirichlet problem is well-posed in 
fl for any 1 < p < oo. In [35], they extended this result to the complements of such 
domains. Their proof used multiple layer potentials introduced by Agmon in [5| in 
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order to solve the Dirichlet problem with continuous boundary data. The general 
outline of their proof parallelled that of the proof of the corresponding result [52] 
for Laplace’s equation. We remark that by [ 1091 Theorem 6.30], we may weaken the 
condition that fl be C 1 to the condition that the unit outward normal v to S2 lies in 
VMO{d£l). (Recall that this condition has been used in [105] ; see formula (13.401) 
above and preceding remarks. This condition was also used in [93]; see Section [572] ) 

As in the case of Laplace’s equation, a result in Lipschitz domains soon followed. 
In [4Tj, Dahlberg, Kenig and Verchota showed that the L p -Dirichlet problem for 
the biharmonic equation is well-posed in any bounded simply connected Lipschitz 
domain fi C R", provided 2 — e < p < 2 + e for some £ > 0 depending on the 
domain Q. 

In 1 1391 . Verchota used the construction of 44, to extend Cohen and Gosselin’s 
results from planar C 1 domains to C 1 domains of arbitrary dimension. Thus, the 
L p -Dirichlet problem for the bilaplacian is well-posed for 1 < p < oo in C 1 domains. 

In ] 140j . Verchota showed that the L p -Dirichlet problem for the polyliarnronic 
operator A m could be solved for 2 — e<p<2 + ein starlike Lipschitz domains 
by induction on the exponent to. He simultaneously proved results for the Ir¬ 
regularity problem in the same range; we will thus delay discussion of his methods 
to Section fOl 

All three of the papers [126], [31] and m constructed biharmonic functions as 
potentials. However, the potentials used differ. m constructed their solutions as 


(4.8) u(X) = 



n— 1 

d 2 n F(X-Y)f(Y)da(Y) + J2 

i =1 



did n F{X -Y)R i 9 (Y) da(Y) 


where Ri are the Riesz transforms. Here F(X) is the fundamental solution to the 
biharmonic equation; thus, u is biharmonic in 1" \ l. As in the case of Laplace’s 
equation, well-posedness of the Dirichlet problem follows from the boundedness 
relation ||AT(Vu)|| iP(an) < C'||/|| Z/P(af2) + C\\g\\ LP{dQ) and from invertibility of the 
mapping (f,g) i-A (u| gn ,d„u) on L p (dfl) x L p (dfl) i-A Wf{dQ) x L p (dfl). 

The multiple layer potential of [34] is an operator of the form 


(4.9) Cf{P) = p.v. f C(P, Q)f{Q) da(Q) 

Jon 

where C(P, Q) is a 3 x 3 matrix of kernels, also composed of derivatives of the 
fundamental solution to the biharmonic equation, and / = (/, f x , f y ) is a “compat¬ 
ible triple” of boundary data, that is, an element of W 1,p (dtt) x L p (dfl) x L p (dfl) 
that satisfies d T f = f x r x + f y r y . Thus, the input is essentially a function and its 
gradient, rather than two functions, and the Riesz transforms are not involved. 

The method of [41] is to compose two potentials. First, the function / € L 2 (9H) 
is mapped to its Poisson extension v. Next, u is taken to be the solution of the 
inhomogeneous equation A u{Y) = (n + 2 Y ■ V)v(Y) with u = 0 on dil. If G(X, Y) 
is the Green’s function for A in H and k Y is the harmonic measure density at V, 
we may write the map /A«as 

(4.10) u(X) = [ G(X,Y)(n + 2Y- V) [ k Y (Q) f{Q) da(Q) dY. 

J O JdSl 

Since (n + 2 Y ■ X)v(Y) is harmonic, u is biharmonic, and so u solves the Dirichlet 
problem. 
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4.2. The L p -Dirichlet problem: the summary of known results on well- 
posedness and ill-posedness. Recall that by [140] , the L p -Dirichlet problem is 
well-posed in Lipscliitz domains provided 2 — s < p < 2 + e. As in the case of 
Laplace’s equation (see m, the range p > 2 — e is sharp. That is, for any p < 2 
and any integers m > 2, n > 2, there exists a bounded Lipschitz domain 52 C R™ 
such that the L p -Dirichlet problem for A m is ill-posed in 52. See m Section 5] for 
the case of the biharmonic operator A 2 , and the proof of Theorem 2.1 in [1201 for 
the polyharmonic operator A m . 

The range p < 2 + e is not sharp and has been studied extensively. Proving or 
disproving well-posedness of the L p -Dirichlet problem for p > 2 in general Lipschitz 
domains has been an open question since [41] . and was formally stated as such in 
[66l Problem 3.2.30]. (Earlier in [29] Question 7], the authors had posed the more 
general question of what classes of boundary data give existence and uniqueness of 
solutions.) 

In ( 1171 Theorem 10.7], Pipher and Vercliota constructed Lipschitz domains 52 
such that the L p -Dirichlet problem for A 2 was ill-posed in 52, for any given p > 
6 (in four dimensions) or any given p > 4 (in five or more dimensions). Their 
counterexamples built on the study of solutions near a singular point, in particular 
upon [95] and [§9]. In Il20l . they provided other counterexamples to show that 
the L p -Dirichlet problem for A m is ill-posed, provided p > 2 (n — 1 )/{n — 3) and 
4 < n < 2 m + 1. They remarked that if n > 2m + 1, then ill-posedness follows from 
the results of [96] provided p > 2m/(m — 1). 

The endpoint result at p = oo is the Agmon-Miranda maximum principle (13.211 
discussed above. We remark that if 2 < po < oo, and the L Po -Dirichlet problem 
is well-posed (or (13.21) holds) then by interpolation, the L p -Dirichlet problem is 
well-posed for any 2 < p < p$. 

We shall adopt the following definition (justified by the discussion above). 

Definition 4.11. Suppose that m > 2 and n > 4. Then is defined to be the 
extended real number that satisfies the following properties. If 2 < p < Pm,m then 
the L p -Dirichlet problem for A m is well-posed in any bounded Lipschitz domain 
52 C R". Conversely, if p > p m ,n, then there exists a bounded Lipschitz domain 
52 C R™ such that the L p -Dirichlet problem for A m is ill-posed in 52. Here, well- 
posedness for 1 < p < oo is meant in the sense of (14.51) . and well-posedness for 
p = oo is meant in the sense of the maximum principle (see (14.251) below). 

As in [41| , we expect the range of solvability for any particular Lipschitz domain 
52 to be 2 — e < p < p rn . n + e for some e depending on the Lipschitz character of 52. 

Let us summarize here the results currently known for p m ,n■ More details will 
follow in Section fOl 

For any m > 2, we have that 

• If n = 2 or n = 3, then the L p -Diriclilet problem for A m is well-posed in 

any Lipschitz domain 52 for any 2 < p < oo. amumi) 

• If 4 < n < 2m + 1, then p m ^ n = 2(n — 1 )/(n — 3). f [129 . 120] .1 

• If n = 2m + 2, then Pm,n = 2m /(to — 1) = 2 (n — 2 )/(n— 4). ( 1130 , 96].) 

• If n > 2 to + 3, then 2 (n — l)/(n — 3) < p m ,n < 2m /(to — 1). ( 129 96].) 

The value of p m . n , for n > 2m + 3, is open. 

In the special case of biharmonic functions (to = 2), more is known. 

• P 2,4 = 6, p 2 ,5 = 4, p 2 ,6 = 4, and p 2 ,7 = 4. ([329! and [130]) 
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• If n > 8, then 

4 

2 H - 7 — < P2,n A 4 

n - A n 


where 

n + 10 + 2y / 2(n 2 — n + 2) 

An = ^ ■ 

( m ) 

• If 12 is a C 1 or convex domain of arbitrary dimension, then the L p -Dirichlet 
problem for A 2 is well-posed in 12 for any 1 < p < oo. ([nnumiin].) 

We comment on the nature of ill-posedness. The counterexamples of m and 
m for p < 2 are failures of uniqueness. That is, those counterexamples are 
nonzero functions u, satisfying A m u = 0 in 12, such that d^u = 0 on <912 for 
0 < k < to — 1, and such that lV(V m-1 it) £ L p (<912). 

Observe that if 12 is bounded andp > 2, then L p (<912) C L 2 (dl2). Because the L 2 - 
Dirichlet problem is well-posed, the failure of wcll-posedness for p > 2 can only be 
a failure of the optimal estimate _/V(V m_1 u) £ L p (dQ). That is, if the L p -Dirichlet 
problem for A m is ill-posed in 12, then for some Whitney array / £ WA T ^ n _ 1 {dVl) 
and some g £ L p {di 2), the unique function u that satisfies A m u = 0 in 12, d a u = f a , 
d™~ 1 u = g and 7V(V m_1 u) £ L 2 (dfl) does not satisfy 2V(V m-1 u) £ L p (dfl). 


4.3. The regularity problem and the L p -Dirichlet problem. In this section 
we elaborate on some of the methods used to prove the Dirichlet well-posedness 
results listed above, as well as their historical context. This naturally brings up a 
consideration of a different boundary value problem, the L 9 -regularity problem for 
higher order operators. 

Recall that for second-order equations the regularity problem corresponds to 
finding a solution with prescribed tangential gradient along the boundary. In anal¬ 
ogy, we say that the L 9 -regularity problem for A m is well-posed in 12 if there exists 
a constant C > 0 such that, whenever / £ W4^(i912), there exists a unique function 
u that satisfies 

A m u = 0 in 12, 

d a ii\ gn = f a for all 0 < \a\ < m — 1, 

\\N(n)\\ Lq{dn) <c J2 H v -/-IIl,( 9 o)- 

\ot\=m —1 

There is an important endpoint formulation at q = 1 for the regularity problem. 
We say that the i? 1 -regularity problem is well-posed if there exists a constant 
C > 0 such that, whenever / lies in the Whitney-Hardy space there 

exists a unique function u that satisfies 

A m u = 0 in 12, 

d a u\ dn = f a for all 0 < |a| < m — 1, 

||iV(V m u )|| L 

qan) — C W^rfaWm^n)- 

| Ct | — 771 — 1 

The space H^d 12) is defined as follows. 
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Definition 4.13. We say that a £ WAf l {dVl) is a IL^(SD)-L 9 atom if a is sup¬ 
ported in a ball B(Q, r) D dfl and if 

H V ^ a «Hw(an) < cr(S(Q,r) naff) 1/9-1 . 

\a\=m—l 

If / € WAj n (dn) and there are H^-L 2 atoms at and constants A& £ C such that 

OO 

V r /a = ^2 AfcV r (a fe ) a for all |a| = m - 1 

k =1 

and such that ^ |Afc| < oo, we say that / £ H^fdtt), with ||/|| ff i being the 
smallest | Afc | among all such representations. 

In 140| . Verchota proved well-posedness of the L 2 -Dirichlet problem and the 
L 2 -regularity problem for the polyharmonic operator A m in any bounded starlike 
Lipschitz domain by simultaneous induction. 

The base case m = 1 is valid in all bounded Lipschitz domains by [3B] and j)2|. 
The inductive step is to show that well-posedness for the Dirichlet problem for A m+1 
follows from well-posedness of the lower-order problems. In particular, solutions 
with d a u = f a may be constructed using the regularity problem for A m , and the 
boundary term d™u = g , missing from the regularity data, may be attained using 
the inhomogeneous Dirichlet problem for A m . On the other hand, it was shown 
that the well-posedness for the regularity problem for A m+1 follows from well- 
posedness of the lower-order problems and from the Dirichlet problem for A m+1 , 
in some sense, by realizing the solution to the regularity problem as an integral of 
the solution to the Dirichlet problem. 

As regards a broader range of p and q , Pipher and Verchota showed in m that 
the L p -Dirichlet and L 9 -regularity problems for A 2 are well-posed in all bounded 
Lipschitz domains Slcl 3 , provided 2 < p < oo and 1 < q < 2. Their method relied 
on duality. Using potentials similar to those of HU, they constructed solutions to 
the L 2 -Dirichlet problem in domains above Lipschitz graphs. The core of their 
proof was the invertibility on L 2 (dfl) of a certain potential operator T. They were 
able to show that the invertibility of its adjoint T* on L 2 (9D) implies that the 
L 2 -regularity problem for A 2 is well-posed. Then, using the atomic decomposition 
of Hardy spaces, they analyzed the H 1 -regularity problem. Applying interpolation 
and duality for T* once again, now in the reverse regularity-to-Dirichlet direction, 
the full range for both regularity and Dirichlet problems was recovered in domains 
above graphs. Localization arguments then completed the argument in bounded 
Lipschitz domains. 

In four or more dimensions, further progress relied on the following theorem of 
Shen. 


Theorem 4.14 ( [150] ). Suppose that Q C l n is a Lipschitz domain. The following 
conditions are equivalent. 

• The L p -Dirichlet problem for L is well-posed, where L is a symmetric el¬ 
liptic system of order 2m with real constant coefficients. 

• There exists some constant C > 0 and some p > 2 such that 

(4.15) (-f iV(V m-1 u) p da) <c(V V(V m-1 u) 2 dcr 

\J B(Q,r)ndn J \J B(Q,2r)ndn 


1/2 








HIGHER-ORDER ELLIPTIC EQUATIONS 


27 


holds whenever u is a solution to the L 2 -Dirichlet problem for L in 57, with 
Vzt = 0 on B(Q , 3r) D 957. 


For the polyharmonic operator A m , this theorem was essentially proven in [12'9| . 
Furthermore, the reverse Holder estimate (14.151) withp = 2(n— l)/(n—3) was shown 
to follow from well-posedness of the L 2 -regularity problem. Thus the L p -Dirichlet 
problem is well-posed in bounded Lipscliitz domains in R" for p = 2 (n — 1 )/(n — 
3). By interpolation, and because reverse Holder estimates have self-improving 
properties, well-posedness in the range 2 < p < 2 (n — l)/(n — 3) + £ for any 
particular Lipschtiz domain follows automatically. 

Using regularity estimates and square-function estimates, Shen was able to fur¬ 
ther improve this range of p. He showed that with p = 2 + 4/(n — A), 0 < A < n, 
the reverse Holder estimate (14.151) is true, provided that 

(4.16) [ |V m_1 zt| 2 < C ( 15 ) A f IV™- 1 ^ 2 

JB(Q,r) nn Jb(Q,r) nn 


holds whenever u is a solution to the L 2 -Dirichlet problem in 57 with N{X m ~ 1 u ) € 
L 2 (di 7) and V fe u| B ^ = 0 for all 0 < k < m — 1. 

It is illuminating to observe that the estimates arising in connection with the 
pointwise bounds on the solutions in arbitrary domains (cf. Section 13.11) and the 
Wiener test (cf. Section ITTTdl) . take essentially the form (14.161) . Thus, Theorem 14. 141 
and its relation to (14.161) provide a direct way to transform results regarding lo¬ 
cal boundary regularity of solutions, obtained via the methods underlined in Sec¬ 
tions I3J1 and IQ into well-posedness of the L p -Diric.hlet problem. 

In particular, consider EH Lemma 5]. If u is a solution to A m u = 0 in B(Q, R) fl 
57, where 57 is a Lipschitz domain, then by EH Lemma 5] there is some constant 
Aq > 0 such that 


(4.17) 


sup M 2 < [ \u(X)\ 2 dX 

B(Q,r) nn b Jb{q,r) nn 


provided that r/R is small enough, that u has zero boundary data on B(Q, i?)fli957, 
and where 57 C R" has dimension n = 2m + 1 or n = 2 m + 2, or where m = 2 and 
n = 7 = 2m + 3. (The bound on dimension comes from the requirement that A m 
be positive with weight F ; see equation (13.461) .') 

It is not difficult to see (cf., e.g., ( 130 ) Theorem 2.6]), that (14.171) implies (14.161) 
for some A > n— 2m+2, and thus implies well-posedness of the L p -Diric.hlet problem 
for a certain range of p. This provides an improvement on the results of [1291 in the 
case m = 2 and n = 6 or n = 7, and in the case m > 3 and n = 2m + 2. Shen has 
stated this improvement in 11301 Theorems 1.4 and 1.5]: the L p -Diric.hlet problem 
for A 2 is well-posed for 2 < p < 4 + e in dimensions n = 6 or n = 7, and the 
L p -Diric.hlet problem for A m is well-posed if 2 < p < 2m/(m — 1) + e in dimension 
n = 2 to + 2. 

The method of weighted integral identities, related to positivity with weight F 
(cf. (13.461) 1. can be further finessed in a particular case of the biharmonic equation. 
m uses this method (extending the ideas from [ST]) to show that if n > 8, then 
(14.161) is valid for solutions to A 2 with A = A„, where 

n + 10 + 2sj2{n 2 — n + 2) 

= 


(4.18) 


7 
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We now return to the L 9 -regularity problem. Recall that in [lT~7j . Pipher and 
Verchota showed that if 2 < p < oo and 1/p+l/q < 1, then the L p -Dirichlet problem 
and the L 9 -regularity problem for A 2 are both well-posed in three-dimensional 
Lipschitz domains. They proved this by showing that, in the special case of a 
domain above a Lipschitz graph, there is duality between the L p -Dirichlet and ir¬ 
regularity problems. Such duality results are common. See [69], j!33j . and [70] for 
duality results in the second-order case; although even in that case, duality is not 
always guaranteed. (See [80).) Many of the known results concerning the regularity 
problem for the polyharmonic operator A m are results relating the L p -Dirichlet 
problem to the iT-regularity problem. 

In 107) . I. Mitrea and M. Mitrea showed that if 1 < p < oo and 1/p + l/q = 
1, and if the L 9 -regularity problem for A 2 and the L p -regularity problem for A 
were both well-posed in a particular bounded Lipschitz domain Cl, then the L p - 
Dirichlet problem for A 2 was also well-posed in Cl. They proved this result (in 
arbitrary dimensions) using layer potentials and a Green representation formula 
for biharmonic equations. Observe that the extra requirement of well-posedness for 
the Laplacian is extremely unfortunate, since in bad domains it essentially restricts 
consideration to p < 2 + e and thus does not shed new light on well-posedness in 
the general class of Lipschitz domains. As will be discussed below, later Kilty and 
Shen established an optimal duality result for biharmonic Dirichlet and regularity 
problems. 

Recall that the formula (14.151) provides a necessary and sufficient condition for 
well-posedness of the L p -Dirichlet problem. In [72], Kilty and Shen provided a 
similar condition for the regularity problem. To be precise, they demonstrated 
that if q > 2 and L is a symmetric elliptic system of order 2m with real constant 
coefficients, then the L 9 -regularity problem for L is well-posed if and only if the 
estimate 


(4.19) 


' B(Q,r)r\Cl 


AT(V" 



B(Q,2r)nn 


A(V’ 


\ 1/2 
1 'u) 2 da J 


holds for all points Q £ dCl , all r > 0 small enough, and all solutions u to the 
L 2 -regularity problem with C7 k u \ b ^q 3 r ) n an = 0 for 0 < k < m — 1. Observe that 
(14.191) is identical to (14.151) with p replaced by q and m — 1 replaced by m. 

As a consequence, well-posedness of the L 9 -regularity problem in Cl for certain 
values of q implies well-posedness of the L p -Dirichlet problem for some values of p. 
Specifically, arguments using interior regularity and fractional integral estimates 
(given in [72j Section 5]) show that (14.191) implies (I4.15[) with 1/p = l/q—l/(n—l). 
But recall from 130] that (14.151) holds if and only if the L p -Dirichlet problem for 
L is well-posed in Cl. Thus, if 2 < q < n — 1, and if the L 9 -regularity problem 
for a symmetric elliptic system is well-posed in a Lipschitz domain Cl, then the 
L p -Dirichlet problem for the same system and domain is also well-posed, provided 
2 < p < p 0 + e where 1/po = l/q — l/(n — 1). 

For the bilaplacian, a full duality result is known. In im Kilty and Shen showed 
that, if 1 < p < oo and l/p+ l/q = 1, then well-posedness of the L p -Dirichlet 
problem for A 2 in a Lipschitz domain Cl, and well-posedness of the L 9 -regularity 
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problem for A 2 in 12, were both equivalent to the bilinear estimate 


(4.20) 


/ 

Jn 


An Av 


< C (||v T v/|| iP + idq-^-^llv/IL. + l^r 2/( "- 1) ||/|| LP ) 
x (nv. 9 |i i9 + i^r i/(ri - i) ii 5 |i i ,) 


for all /, g £ C'o°(]R n ), where u and v are solutions of the L 2 -regularity problem 
with boundary data d a u = d a f and d a v = d a g. Thus, if 12 C R" is a bounded 
Lipschitz domain, and if l/p + 1/q = 1, then the Z/'-Dirichlet problem is well-posed 
in 12 if and only if the L 9 -regularity problem is well-posed in 12. 

All in all, we see that the L p -regularity problem for A 2 is well-posed in 12 C R” 
if 

• 12 is C 1 or convex, and 1 < p < oo. 

• n = 2 or n = 3 and 1 < p < 2 + e. 

• n = 4 and 6/5 — e < p < 2 + e. 

• n = 5, 6 or 7, and 4/3 — e < p <2 + e. 

• n > 8, and 2 — 4+ra 4 _ A <p < 2 + e, where \ n is given by (14.181) . The above 
ranges of p are sharp, but this range is still open. 


4.4. Higher-order elliptic systems. The polyharmonic operator A m is part of 
a larger class of elliptic higher-order operators. Some study has been made of 
boundary-value problems for such operators and systems. 

The L p -Dirichlet problem for a strongly elliptic system L of order 2m, as defined 
in Definition 12.11 is well-posed in 12 if there exists a constant C such that, for 
every / £ Wb4(( l _ 1 (912 i-a C/) and every g £ L p (<912 i-a C £ ), there exists a unique 
vector-valued function u : 12 <—>■ C e such that 


(4.21) 


( L u)j = 0 

k —1 \a\=\/3\=m 

9 a U = f a 
d^u = g 


in 12 for each 1 < j < t, 

on <912 for |a| < m — 2, 
on <912, 


l|JV(V”>- 1 »)||„ (M) <c Y. 11 Vr/a 11^5(90) +c\m Lp(dny 

\a\=m —2 


The L 9 -regularity problem is well-posed in 12 if there is some constant C such that, 
for every / £ WA p t (<912 i-a C £ ), there exists a unique u such that 


(4.22) 


{Lu)j = ^ ^ d a a^ l3 d l3 Uk = 0 in 12 for each 1 < j < £, 

k= 1 \a\=\/3\=m 

d a u = f a on <912 for |a| < m — 1, 
\\N(V m u)\\ L , {dn) <C £ 


|o!|=m—1 


In m, Pipher and Verchota showed that the L p -Dirichlet and L p -regularity 
problems were well-posed for 2 — e < p < 2 + e, for any higher-order elliptic 
partial differential equation with real constant coefficients, in Lipschitz domains of 
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arbitrary dimension. This was extended to symmetric elliptic systems in [141] . A 
key ingredient of the proof was the boundary Garding inequality 

7 [ \X m u\(-v n )da 
4 Jd n 

£ 

< T, V [ d a a^gd 0 u k {-i'n) da + C f \X7 m ~ 1 d n ufda 

j,k=l |a| = |/3|=rn 9f2 ■' dQ 

valid if u G Cg°(R n )*> if L = d a a%d^ is a symmetric elliptic system with real 
constant coefficients, and if 12 is the domain above the graph of a Lipschitz function. 
We observe that in this case, (— v n ) is a positive number bounded from below. 
Pipher and Verchota then used this Garding inequality and a Green’s formula to 
construct the nontangential maximal estimate. See [120 : and }f4ll Sections 4 and 6 ]. 

As in the case of the poly harmonic operator A m , this first result concerned the 
P p -Dirichlet problem and L 9 -regularity problem only for 2 — e < p < 2 + £ and for 
2 — s < q < 2 + e. The polyharmonic operator A m is an elliptic system, and so we 
cannot in general improve upon the requirement that 2 — e < p for well-posedness 
of the L p -Dirichlet problem. 

However, we can improve on the requirement p < 2+e. Recall that Theorem l4.14l 
from |130j , and its equivalence to (14.161) , were proven in the general case of strongly 
elliptic systems with real symmetric constant coefficients. As in the case of the 
polyharmonic operator A m , (14.1511 follows from well-posedness of the P 2 -regularity 
problem provided p = 2(n—l)/(n—3), and so if L is such a system, the L p -Dirichlet 
problem for L is well-posed in 11 provided 2 — e < p < 2(n — l)/(n — 3) + £. This is 
[Ml Corollary 1.3]. Again, by the counterexamples of [ 120] . this range cannot be 
improved if m > 2 and 4 < n < 2m + 1; the question of whether this range can be 
improved for general operators L if n > 2m + 2 is still open. 

Little is known concerning the regularity problem in a broader range of p. Recall 
that (14.191) from [72] was proven in the general case of strongly elliptic systems with 
real symmetric constant coefficients. Thus, we known that for such systems, well- 
posedness of the L 9 -regularity problem for 2 < q < n — 1 implies well-posedness of 
the L p -Dirichlet problem for appropriate p. The question of whether the reverse 
implication holds, or whether this result can be extended to a broader range of q , 
is open. 

4.5. The area integral. One of major tools in the theory of second-order elliptic 
differential equations is the Lusin area integral, defined as follows. If w lies in 
Wi 2 j oc (H) for some domain H C R ra , then the area integral (or square function) of 
w is defined for Q G dtt as 

Sw(Q)=\ [ \S7w(X)\ 2 dist{X,dn) 2 ~ n dX 

\Jr(Q) 

In [37) . Dahlberg showed that if u is harmonic in a bounded Lipschitz domain 12, 
if Po G 12 and u(Pq) = 0, then for any 0 < p < oo, 

(4.23) i f Su p da < f ( Nu) p da < C [ ( Su) p da 

k Jon Jon Jon 

for some constants C depending only on p , 12 and Pq. Thus, the Lusin area integral 
bears deep connections to the P p -Dirichlet problem. In [38], Dahlberg, Jerison and 
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Kenig generalized this result to solutions to second-order divergence-form elliptic 
equations with real coefficients for which the ZA-Dirichlet problem is well-posed for 
at least one r. 

If L is an operator of order 2 to, then the appropriate estimate is 

(4.24) I/" N{Vu m ~ 1 ) p da < f S(Vu rn ~ 1 ) p da < C f N{Vu m ~ 1 ) p da. 

^ Jon Jd n Jo n 

Before discussing their validity for particular operators, let us point out that such 
square-function estimates are very useful in the study of higher-order equations. 
In 1 :T0 . Shen used (14.241) to prove the equivalence of (14.161) and (14.151) . above. In 
d], Kilty and Shen used (14.241) to prove that well-posedness of the L p -Dirichlet 
problem for A 2 implies the bilinear estimate (14.201) . The proof of the maximum 
principle (13.21) in 1141 , Section 8 ] (to be discussed in Section 14.61) also exploited 
(14.241) . Estimates on square functions can be used to derive estimates on Besov 
space norms; see [H Proposition S]. 

In 1 116j . Pipher and Verchota proved that (14.241) (with m = 2) holds for solutions 
u to A 2 it = 0, provided is a bounded Lipschitz domain, 0 < p < oo, and 
Vu(-Po) = 0 for some fixed Pq £ 11. Their proof was an adaptation of Dahlberg’s 
proof [37j of the corresponding result for harmonic functions. They used the L 2 - 
tlreory for the biharmonic operator EH, the representation formula (14.101) , and the 
L 2 -theory for harmonic functions to prove good-A inequalities, which, in turn, imply 
L p estimates for 0 < p < oo. 

In [40], Dahlberg, Kenig, Pipher and Verchota proved that (14.241) held for solu¬ 
tions u to Lu = 0, for a symmetric elliptic system L of order 2 m with real constant 
coefficients, provided as usual that 11 is a bounded Lipschitz domain, 0 < p < oo, 
and V m_ 1 u(Po) = 0 for some fixed Pq £ H. The argument is necessarily consid¬ 
erably more involved than the argument of II 16| or m- In particular, the bound 
|| 5 '(V ,n_ 1 u)|| i 2 ( a o) < C'||lV(V ,n_ 1 M)Hz j 2 (£) 0 ) was proven in three steps. 

The first step was to reduce from the elliptic system L of order 2m to the scalar 
elliptic operator M = det L of order 2£m, where i is as in formula (12.21) . The 


second step was to reduce to elliptic equations of the form Y^\ a \=m a a d 2a u = 0, 
where |a a | > 0 for all |a| = m. Finally, it was shown that for operators of this form 


V [ a a d a u(X) 2 dist(A,5H)dA < C [ iV(V Tn_1 «) 
.“1 J n Jan 


da. 


\ot\=m 


The passage to 0 < p < oo in (14.241) was done, as usual, using good-A inequalities. 
We remark that these arguments used the result of 119 that the L 2 -Diric.hlet 
problem is well-posed for such operators L in Lipschitz domains. 

It is quite interesting that for second-order elliptic systems, the only currently 
known approach to the square-function estimate (14.231) is this reduction to a higher- 
order operator. 


4.6. The maximum principle in Lipschitz domains. We are now in a posi¬ 
tion to discuss the maximum principle (13.21) for higher-order equations in Lipschitz 
domains. 

We say that the maximum principle for an operator L of order 2m holds in the 
bounded Lipschitz domain Q if there exists a constant C > 0 such that, whenever 
/ £ WA'^_ l (d^l) C WA^ n _ l {dVL) and g £ L°°(dfl) C L 2 (dfl), the solution u to the 
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Dirichlet problem (14.211) with boundary data / and g satisfies 

(4.25) ||V m - 1 W || ioo <C'|| ff || iao(an) +C £ ||V r /J Loo(sn) . 

\a\=m—2 

The maximum principle (14.251) was proven to hold in three-dimensional Lips- 
chitz domains bv Pipher and Verchota in I118| (for biharmonic functions), in m 
(for polyharmonic functions), and by Verchota in (1411 Section 8] (for solutions 
to symmetric systems with real constant coefficients). Pipher and Verchota also 
proved in 11181 that the maximum principle was valid for biharmonic functions in 
C 1 domains of arbitrary dimension. In m Theorem 1.5], Kilty and Shen observed 
that the same techinque gives validity of the maximum principle for biharmonic 
functions in convex domains of arbitrary dimension. 

The proof of |118j uses the L 2 -regularity problem in the domain 12 to construct 
the Green’s function G(X,Y) for A 2 in 12. Then if u is biharmonic in 12 with 
IV(Vzt) £ L 2 (d 12), we have that 

(4.26) u(X) = [ u(Q) d„AG(X, Q) da(Q) + [ d v u{Q) AG(A,Q) da(Q) 

J dCt J dQ, 

where all derivatives of G are taken in the second variable Q. If the H 1 -regularity 
problem is well-posed in appropriate subdomains of 12, then V 2 VxG(A', •) is in 
L 1 (912) with i 1 -norm independent of A', and so the second integral is at most 
C'||<9 ! /'u|| Z/00 ( an ). By taking Riesz transforms, the normal derivative d„AG(X,Q) 
may be transformed to tangential derivatives V r AG(A, Q); integrating by parts 
transfers these derivatives to u. The square-function estimate (14.241) implies that 
the Riesz transforms of X xAqG(X, Q) are bounded on L 1 (i912). This completes 
the proof of the maximum principle. 

Similar arguments show that the maximum principle is valid for more general 
operators. See [120] for the polyharmonic operator, or [1411 Section 8] for arbitrary 
symmmetric operators with real constant coefficients. 

An important transitional step is the wcll-posedness of the H 1 -regularity prob¬ 
lem. It was established in three-dimensional (or C 1 ) domains in [118i Theorem 4.2] 
and 11201 Theorem 1.2] and discussed in [1411 Section 7]. In each case, well- 
posedness was proven by analyzing solutions with atomic data / using a technique 
from [39]. A crucial ingredient in this technique is the well-posedness of the L p - 
Dirichlet problem for some p < (n — l)/(n — 2); the latter is valid if n = 3 by 
[44]. and (for A 2 ) in G 1 and convex domains bv [140] and EH, but fails in general 
Lipschitz domains for n > 4. 

4.7. Biharmonic functions in convex domains. We say that a domain 12 is 
convex if, whenever A, Y £ 12, the line segment connecting X and Y lies in 12. 
Observe that all convex domains are necessarily Lipschitz domains but the converse 
does not hold. Moreover, while convex domains are in general no smoother than 
Lipschitz domains, the extra geometrical structure often allows for considerably 
stronger results. 

Recall that in [82j, the second author of this paper and Maz’ya showed that 
the gradient of a biharmonic function is bounded in a three-dimensional domain. 
This is a sharp property in dimension three, and in higher dimensional domains the 
solutions can be even less regular (cf. Section I3.1D . However, using some intricate 
linear combination of weighted integrals, the same authors showed in [81] that 
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second derivatives to biharmonic functions were locally bounded when the domain 
was convex. To be precise, they showed that if 11 is convex, and u G W%{£1) is a 
solution to A 2 u = h for some h G Cg°(fl \ B{Q , 101?)), R > 0, Q G dfl, then 

(4.27) sup |V 2 u|<-^(-/ 

B(Q,R/ 5)nn n \Jnr\B{Q,5R)\B(Q,R/2) 

In particular, not only are all boundary points of convex domains 1-regular, but 
the gradient Vu is Lipschitz continuous near such points. 

Kilty and Shen noted in m that (14.2711 implies that (14.1911 holds in convex 
domains for any q; thus, the L 9 -regularity problem for the bilaplacian is well-posed 
for any 2 < q < oo in a convex domain. Well-posedness of the L p -Dirichlet problem 
for 2 < p < oo has been established by Shen in [131] . By the duality result (14.201) , 
again from [ZD, this implies that both the L p -Dirichlet and L 9 -regularity problems 
are well-posed, for any 1 < p < oo and any 1 < q < oo, in a convex domain of 
arbitrary dimension. They also observed that, by the techniques of (1181 (discussed 
in Section 14.61 above), the maximum principle (14.251) is valid in arbitrary convex 
domains. 

It is interesting to note how, once again, the methods and results related to 
pointwise estimates, the Wiener criterion, and local regularity estimates near the 
boundary are intertwined with the well-posedness of boundary problems in L p . 

4.8. The Neumann problem for the biharmonic equation. So far we have 
only discussed the Dirichlet and regularity problems for higher order operators. 
Another common and important boundary-value problem that arises in applica¬ 
tions is the Neumann problem. Indeed, the principal physical motivation for the 
inhomogeneous biharmonic equation A 2 u = h is that it describes the equilibrium 
position of a thin elastic plate subject to a vertical force h. The Dirichlet problem 
u \dn = Vu| an = g describes an elastic plate whose edges are clamped, that is, 
held at a fixed position in a fixed orientation. The Neumann problem, on the other 
hand, corresponds to the case of a free boundary. Guido Sweers has written an 
excellent short paper 1371 discussing the boundary conditions that correspond to 
these and other physical situations. 

More precisely, if a thin two-dimensional plate is subject to a force h and the 
edges are free to move, then its displacement u satisfies the boundary value problem 

A 2 u = h in fi, 
pAu + (1 — p)d1u = 0 on c>f2, 

<9„A u + (1 — p)d TTV u = 0 on dfl. 

Here p is a physical constant, called the Poisson ratio. This formulation goes back 
to Kirchoff and is well known in the theory of elasticity; see, for example, Section 3.1 
and Chapter 8 of the classic engineering text m- We remark that by i 1 13. formula 
(8-10)], 

d„Au + (1 - p)d TTV u = d„Au + (1 - p)d T {d VT u). 

This suggests the following homogeneous boundary value problem in a Lipschitz 
domain H of arbitrary dimension. We say that the L p -Neumann problem is well- 
posed if there exists a constant C > 0 such that, for every f 0 € L p (dfl) and 
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Aq € Wf_ 1 (dCl), there exists a function u such that 


(4.28) 


A 2 u = 0 

M p u := pAu + (1 - p)dlu = f 0 


in fl, 
on dfl. 


K p u := d„Au + (1 - p)-d Tij {d VTij u) = A 0 on <90, 
_ l|A r (V 2 u)|| LP(d Q) — C|l/o|liYf(an) + C||Ao|lw^ i(an)- 


Here = fje-j — Vje.i is a vector orthogonal to the outward normal v and lying in 
the XiXj-p\ane. 

In addition to the connection to the theory of elasticity, this problem is of interest 
because it is in some sense adjoint to the Dirichlet problem (14.41) . That is, if 
A 2 u = A 2 w = 0 in O, then f gQ d„w M p u — w K p u da = J gfl d v u M p w — u K p w da, 
where M p and K p are as in (14.281) : this follows from the more general formula 


(4.29) / wA 2 u= / (pAu Aw + (l — p)djkudjkw) + / < 

JQ, JQ, J d£l 


w KpU — dviu MpU da 


valid for arbitrary smooth functions. 
Green’s identity for the Laplacian 


This formula is analogous to the classical 


(4.30) 


/ 

Jn 


> Art = — / Vu-Vw 


/ 

J a n 


w v ■ Vm da. 


Observe that, contrary to the Laplacian or more general second order operators, 
there is a family of relevant Neumann data for the biharmonic equation. Moreover, 
different values (or, rather, ranges) of p correspond to different natural physical 
situations. We refer the reader to 1142 for a detailed discussion. 

In 2^, Cohen and Gosselin showed that the L p -Neumann problem (14.281) was 
well-posed in C ' 1 domains contained in K 2 for for 1 < p < oo, provided in addition 
that p = — 1. The method of proof was as follows. Recall from (14.91) that Cohen 
and Gosselin showed that the L p -Dirichlet problem was well-posed by constructing 
a multiple layer potential Cf with boundary values (/ + /C)/, and showing that 
I + 1C is invertible. We remark that because Cohen and Gosselin preferred to work 
with Dirichlet boundary data of the form [u,d x u,d y u)\ g ^ rather than of the form 
(u,d v u) the notation of [55| is somewhat different from that of the present 
paper. In the notation of the present paper, the method of proof of [35) was to 
observe that (/ + K,)*9 is equivalent to {K_iv6, M_iv0)qqc , where v is another 
biharmonic layer potential and (/ + 1C)* is the adjoint to (/ + JC). Well-posedness 
of the Neumann problem then follows from invertibility of I + K. on dCL c . 

In |142j . Verchota investigated the Neumann problem (14.281) in full generality. 
He considered Lipschitz domains with compact, connected boundary contained in 
R", n > 2. He showed that if — l/(n — 1) < p < 1, then the Neumann problem is 
well-posed provided 2 —e < p < 2+e. That is, the solutions exist, satisfy the desired 
estimates, and are unique either modulo functions of an appropriate class, or (in 
the case where H is unbounded) when subject to an appropriate growth condition. 
See ll42l Theorems 13.2 and 15.4]. Verchota’s proof also used boundedness and 
invertibility of certain potentials on L p (dfl); a crucial step was a coercivity estimate 
||V 2 Lt|| i 2 ( £)0 ) < C\\K p u\\ W 2 ^ gn) + C\\M p u\\ L 2 ^ gn y (This estimate is valid provided 
u is biharmonic and satisfies some mean-value hypotheses; see [142 . Theorem 7.6]). 
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More recently, in [1 132 . , Shen improved upon Verchota’s results by extending the 
range on p (in bounded simply connected Lipschitz domains) to 2(n— l)/(n+l)—£ < 
p<2 + eifn>4, and l<p<2 + eifn = 2orn = 3. This result again was 
proven by inverting layer potentials. Observe that the L p -regularity problem is also 
known to be well-posed for p in this range, and (if n > 6) in a broader range of p; 
see Section 14.31 The question of the sharp range of p for which the L p -Neumann 
problem is well-posed in a Lipschitz domain is still open. 

Finally, in [1091 Section 6.5], I. Mitrea and M. Mitrea showed that if Q C K" 
is a simply connected domain whose unit outward normal v lies in VMO(dfl) (for 
example, if fl is a C 1 domain), then the acceptable range of p is 1 < p < oo; 
this may be seen as a generalization of the result of Cohen and Gosselin to higher 
dimensions, to other values of p, and to slightly rougher domains. 

It turns out that extending the well-posedness results for the Neumann problem 
beyond the case of the bilaplacian is an excruciatingly difficult problem, even if 
one considers only fourth-order operators with constant coefficients. Even defining 
Neumann boundary values for more general operators is a difficult problem (see 
Section m , and while some progress has been made (see [?,, 109 - 22 20], or Sec¬ 
tion [53] below), at present there are no well-posedness results for the Neumann 
problem with L p boundary data. 

In analogy to (14.291) and (14.301) , one can write 

(4.31) / w Lu = A[u , w]+ w Kau — d v w Mau da, 

Jn Jd n 

where A[u,w\ = J]|a|=|/ 3|=2 °a /3 is an energy form associated to the 

operator L = ]C| a |=|/3|=2 a a pD a D 13 . Note that in the context of fourth-order op¬ 
erators, the pair {w, d u w) constitutes the Dirichlet data for w on the boundary, 
and so one can say that the operators Kau and Mau define the Neumann data 
for u. One immediately faces the problem that the same higher-order operator L 
can be rewritten in many different ways and gives rise to different energy forms. 
The corresponding Neumann data will be different. (This is the reason why there 
is a family of Neumann data for the biharmonic operator.) 

Furthermore, whatever the choice of the form, in order to establish well-posedness 
of the Neumann problem, one needs to be able to estimate all second derivatives of 
a solution on the boundary in terms of the Neumann data. In the analogous second- 
order case, such an estimate is provided by the Rellich identity, which shows that 
the tangential derivatives are equivalent to the normal derivative in L 2 for solu¬ 
tions of elliptic PDEs. In the higher-order scenario, such a result calls for certain 
coercivity estimates which are still rather poorly understood. We refer the reader 
to [143 for a detailed discussion of related results and problems. 


4.9. Inhomogeneous problems and other classes of boundary data. In [4], 
Adolfsson and Pipher investigated the inhomogeneous Dirichlet problem for the 
biharmonic equation with data in Besov and Sobolev spaces. While resting on the 
results for homogeneous boundary value problems discussed in Sections 14. ll ancl Td. 31 
such a framework presents a completely new setting, allowing for the inhomogeneous 
problem and for consideration of classes of boundary data which are, in some sense, 
intermediate between the Dirichlet and the regularity problems. 
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They showed that if / G WA\ +s (dVL) and h £ Z/ +]/p _ 3 (S4), then there exists a 
unique function u that satisfies 


(4.32) 


| A 2 u = h in fi, 

( Tr d a u = f a , for 0 < |a| < 1 

subject to the estimate 

(4.33) IMIi* 1 /;p+ 1 (n) - 17 («) + C \\f\\wA p 1+s (dQ) 


provided 2 — e < p < 2 + e and 0 < s < 1. Here Tr w denotes the trace of w in the 
sense of Sobolev spaces; that these may be extended to functions u £ L^ +1+1 ^ p , 
s > 0, was proven in [4] Theorem 1.12]. 

In Lipschitz domains contained in R , they proved these results for a broader 
range of p and s, namely for 0 < s < 1 and for 


(4.34) 


max 1 , 


s + 1 + £ 


< p < 



S < £, 

£ < S < 1 . 


Finally, in C 1 domains, they proved these results for any p and s with 1 < p < oo 
and 0 < s < 1. 

In 110 . I. Mitrea, M. Mitrea and Wright extended the three-dimensional results 
to p = oo (for 0 < s < e) or 2/(s + l+£)<p<l (for 1 — e < s < 1). They also 
extended these results to data h and / in more general Besov or Triebel-Lizorkin 
spaces. 

In )108| . I. Mitrea and M. Mitrea extended the results of [4] to higher dimensions. 
That is, they showed that if f2 C R" is a Lipschitz domain and n > 4, then there 
is a unique solution to the problem (14.321) subject to the estimate (14.331) provided 
that 0 < s < 1 and that 


max 


n — 1 


s + (n - l )/2 + £ 


< p < 


00 , 


n— 1 


(n—3)/2+s—£ ’ 


(n — 3)/2 + s < £, 
£ < s < 1 . 


As in m, their results extend to more general function spaces. 

I. Mitrea and M. Mitrea also showed that, for the same values of p and s, there 
exist unique solutions to the inhomogeneous Neumann problem 


{ A 2 u = h in H, 

M p u = f on dfl, 

K p u = A on dfl 

where M p and K p are as in Section I4~51 subject to the estimate 


(4.35) H w llLj + 1 /p+ 1 (f2) — ^1I^II.L* + 1 /p _3(n) + C'll/llsfy’^an) + 


Finally, in 11091 Section 6.4], I. Mitrea and M. Mitrea proved similar results for 
more general constant-coefficient elliptic operators. That is, if L is an operator 
given by formula (12.21) whose coefficients satisfy the ellipticity condition (12.31) . and 
if fl C R™ is a bounded Lipschitz domain, then there exists a unique solution to 
the Dirichlet problem 



in fl, 

for 0 < M < 1 


(4.36) 
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subject to the estimate 




(4.37) IMIs^ , . (a) - Cll^ll B p ’ q (fi) + C\\f\\ W A p ’ q 

provided 2 — e < p < 2 + s, 2 — e < q < 2 + e, and 1/2 — e < s < 1/2 + e. (A very 
similar result is valid for variable-coefficient divergence-form operators; see m, 
discussed as formula (15.91) below.) Furthermore, with a slightly stronger ellipticity 
condition 


Re 


£ £ > aici= 

\a\ — \P\=m j,k=l 


they established well-posedness of the inhomogeneous Neumann problem for ar¬ 
bitrary constant-coefficient operators. See Section 15.51 below for a formulation of 
Neumann boundary data in the case of arbitrary operators. Finally, if L is self- 
adjoint and n > 2m, and if the unit outward normal v to dfl lies in VMO(dfl), then 
the Dirichlet problem (14.361) has a unique solution satisfying the estimate (14.371) for 
any 0 < s < 1 and any 1 < p < oo, 1 < q < oo. 

Let us define the function spaces appearing above. L£(K") is defined to be 
{g : (I — A) a l 2 g £ L P (K”)}; we say g £ (D) if g = h\ n for some h £ L(((R ra ). If k 

is a nonnegative integer, then L p k = W p . If m is an integer and 0 < s < 1, then the 
Whitney-Besov space WA p m _ 1+s = WA p r ] p _ 1+s or WA^_ l+s is defined analogously 
to WA p m (see Definition 14. 6p . except that we take the completion with respect to 
the Whitney-Besov norm 


(4.38) £ ll^ll„ (an) + £ iioiIb.- q (dQ) 

\a\<m—l |cn|—m—1 

rather than the Whitney-Sobolev norm 

E ll^ll L - ( an)+ E H V ^II LP(dtt)- 

\a\<m— 1 \a\=m—l 


In [ij, the general problem (14.321) for A 2 was first reduced to the case h = 0 
(that is, to a homogeneous problem) by means of trace/extension theorems, that 
is, subtracting w(X) = f Rn F(X, Y) h(Y ) dY, and showing that if h £ 7/ + 1 / , p _ 3 (fI) 
then (TruqTrVw;) £ WA P +S (dtt). Next, the well-posedness of Dirichlet and regu¬ 
larity problems discussed in Sections 14.11 and 14.31 provide the endpoint cases s = 0 
and s = 1, respectively. The core of the matter is to show that, if u is biliarnronic, 
k is an integer and 0 < a < 1, then u £ L k+a (il) if and only if 


(4.39) f |V fc+ 1 u(A)| P dist(A', dfl) p ~ pa + |V fe u(A)| P + \u(X)\ p dX < oo, 

■In 

(cf. [U Proposition S]). With this at hand, one can use square-function estimates 
to justify the aforementioned endpoint results. Indeed, observe that for p = 2 the 
first integral on the left-hand side of (14.391) is exactly the L 2 norm of S'(V fe u). 
The latter, by 1 lii (discussed in Section ll.Sp , is equivalent to the L 2 norm of 
the corresponding non-tangential maximal function, connecting the estimate (14.331) 
to the nontangential estimates in the Dirichlet problem (14.41) and the regularity 
problem 14.121 Finally, one can build an interpolation-type scheme to pass to well- 
posedness in intermediate Besov and Sobolev spaces. 

The solution in il09i to the problem (14.361) . at least in the case of general Lip- 
schitz domains, was constructed in the opposite way, by first reducing to the case 
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where the boundary data / = 0. Using duality it is straightforward to establish 
well-posedness in the case p = q = 2 , s = 1 / 2 ; perturbative results then suffice to 
extend to p , q near 2 and s near 1 / 2 . 


5. Boundary value problems with variable coefficients 


Results for higher order differential equations with variable coefficients are very 
scarce. As we discussed in Section^ there are two natural manifestations of higher- 
order operators with variable coefficients. Operators in divergence form arise via the 
weak formulation framework. Conversely, operators in composition form generalize 
the bilaplacian under a pull-back of a Lipschitz domain to the upper half-space. 

Both classes of operators have been investigated. However, operators in diver¬ 
gence form have received somewhat more study; thus, we begin this section by 
reviewing the definition of divergence-form operator. A divergence-form operator 
A, acting on ioc (f2 ^ C £ ), may be defined weakly via (Ell); we say that Lu = h 
if 


(5.1) 




d p u k 


for all tp smooth and compactly supported in f 1 . 


5.1. The Kato problem and the Riesz transforms. We begin with the Kato 
problem and the properties of the Riesz transform; this is an important topic in 
elliptic theory, which formally stands somewhat apart from the well-posedness is¬ 
sues. 

Suppose that L is a variable-coefficient operator in divergence form, that is, an 
operator defined by EH). Suppose that L satisfies the bound 


(5.2) 


H [ a % dP h d<x 9j 

\a\=\P\=mj,k=l JRn 


— C'||V m /|| i 2( R n)||V m (l|| i 2 (R n ) 


for all / and g in W^(M n ), and the ellipticity estimate 



for all functions tp £ Co°(f2 K>• C/). Notice that this is a weaker requirement than 
the pointwise ellipticity condition (12.61) . 

Auscher, Hofmann, McIntosh and Tchamitchian m proved that under these 
conditions, the Kato estimate 

(5-4) ^l|V m /|| L 2 (K „) < \\VLf\\ L2{Rn) < C||V m /|| i2(R „ ) 

is valid for some constant C. They also proved similar results for operators with 
lower-order terms. 

It was later observed in El that by the methods of [18], if 1 < n < 2m, then 
the bound on the Riesz transform V m A -1 / 2 in L p (that is, the first inequality in 
formula (15.41) 1 extends to the range 1 < p < 2 + e, and the reverse Riesz transform 
bound (that is, the second inequality in formula (15.41) 1 extends to the range 1 < p < 
oo. This also holds if the Schwartz kernel Wt.(X,Y) of the operator e~ tL satisfies 
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certain pointwise bounds (e.g., if L is second-order and the coefficients of A are 
real). 

In the case where n > 2m, the inequality ||V m L^ 1 ' /2 /|| iP Q R „) < C’||/Hlp(r«) 
holds for 2n/(n + 2m) — e < p < 2; see [26] ]TT]. The reverse inequality holds 
for max(2 n/(n + 4to) — e, 1) < p < 2 by II, Theorem 18], and for 2 < p < 
2n/(n — 2m) + e by duality (see [12] Section 7.2]). 

In the case of second-order operators, the Kato estimate implies well-posedness 
of boundary-value problems with L 2 data in the upper half space for certain coeffi¬ 
cients in a special (“block”) form. We conjecture that the same is true in the case 
of higher-order operators; see Section EU 


5.2. The Dirichlet problem for operators in divergence form. In this section 
we discuss boundary-value problems for divergence-form operators with variable 
coefficients. At the moment, well-posedness results for such operators are restricted 
in that the boundary problems treated fall strictly between the range of L p -Dirichlet 
and L p -regularity, in the sense of Section 14.91 That is, there are at present no 
well-posedness results for the L p -Dirichlet, regularity, and Neumann problems on 
Lipschitz domains with the usual sharp estimates in terms of the non-tangential 
maximal function for these divergence-form operators. (Such problems are now 
being considered; see Section 15.61 ) 

To be more precise, recall from the discussion in Section 14.91 that the classical 
Dirichlet and regularity problems, with boundary data in L p , can be viewed as the 
s = 0, 1 endpoints of the boundary problem studied in [4], [110 1 and 1108] 

A 2 it = h in Cl, d a u\ gn = f a for all |a| < 1 

with / lying in an intermediate smoothness space WA\ +s (dCl), 0 < s < 1. In the 
context of divergence-form higher-order operators with variable coefficients, essen¬ 
tially the known results pertain only to boundary data of intermediate smoothness. 

In [7], Agranovich investigated the inhomogeneous Dirichlet problem, in Lips¬ 
chitz domains, for such operators L that are elliptic (in the pointwise sense of (12.61) . 
and not the more general condition (15.31) ) and whose coefficients are Lipschitz 
continuous in Cl. 

He showed that if h G L P _ m _ 1+1 i p+s {Cl) and / S WA p m _ 1+s (dCl ), for some 0 < 
s < 1, and if \p — 2| is small enough, then the Dirichlet problem 


(5.5) 



in Cl, 

for all 0 < |a| < m — 1 


has a unique solution u that satisfies the estimate 


(5.6) 


\\ u \\L p 


,-l + s + l/p 


(ft) S 






tS' / 


Agranovich also considered the Neumann problem for such operators. As we 
discussed in Section 14.81 defining the Neumann problem is a delicate matter. In 
the context of zero boundary data, the situation is a little simpler as one can take 
a formal functional analytic point of view and avoid to some extent the discussion 
of estimates at the boundary. We say that u solves the Neumann problem for L, 
with homogeneous boundary data, if the equation EH) in the weak formulation of 
L is valid for all test functions p compactly supported in R” (but not necessarily 
in Cl.) Agranovich showed that, if h € L p _ m _ 1+1 / p+s (Cl), then there exists a unique 
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function u £ L P m _ l+x / p+g (f2) that solves this Neumann problem with homogeneous 
boundary data, under the same conditions on p, s, L as for his results for the 
Dirichlet problem. Here h £ L^(Cl) if h = g | n for some g £ L£,(R") that in 
addition is supported in Cl. 

In 02, Maz’ya, M. Mitrea and Slraposhnikova considered the Dirichlet problem, 
again with boundary data in intermediate Besov spaces, for much rougher coeffi¬ 
cients. They showed that if / £ WAf n _ l+a {dCl), for some 0 < s < 1 and some 
1 < p < oo, if h lies in an appropriate space, and if L is a divergence-form operator 
of order 2m (as defined by (15.11) 1. then under some conditions, there is a unique 
function u that satisfies the Dirichlet problem (15.51) subject to the estimate 


(5.7) 


ir 


,1 — s — l/p 


( [ |d“u(X)| p dist(X,<9D) p - ps - 1 

da|<-m n 


\ !/P 

dX | < OO. 


See [93] Theorem 8.1]. The inhomogeneous data h is required to lie in the space 

the dual S P &Ce t0 KL+1/P— Where 


(5.8) |Mll 


|cx| <m 1 


\ i /p 

\d a u{X)\ p dist(X, dn)P a +P^-P m dX 1 . 


Notice that w £ a if and only if w £ a and d a w = 0 on dCl for all 0 < |a| < 
to — 1 . 

The conditions are that the coefficients a satisfy the weak ellipticity condition 
(15.31) considered in the theory of the Kato problem, that fl be a Lipschitz domain 
whose normal vector v lies in VMO(dfl), and that the coefficients lie in L°°(R n ) 
and in VMO(R n ). Recall that this condition on fi has also arisen in [1051 (it ensures 
the validity of formula (13.401) ). Notice that the L°° bound on the coefficients is 
a stronger condition than the bound 15.21 of US], and the requirement that the 
coefficients lie in VMO(M. n ) is a regularity requirement that is weaker than the 
requirement of [7] that the coefficients be Lipschitz continuous. 

In fact, 93] provides a more intricate result, allowing one to deduce a well- 
posedness range of s and p, given information about the oscillation of the coeffi¬ 
cients and the normal to the domain v. In the extreme case, when the oscilla¬ 
tions for both are vanishing, the allowable range expands toO<s<l,l<p<oo, 
as stated above. 

The construction of solutions to the Dirichlet problem may be simplified using 
trace and extension theorems. In [S3] Proposition 7.3], the authors showed that if 
/ £ WA p m _ 1+s {dVL), then there exists a function F £ a such that d a F = f a on 
dtt. It is easy to see that if F £ W^ a , and the coefficients of L are bounded 
pointwise, then LF £ R_ m! i_ s _i/ p . Thus, the Dirichlet problem 


Lu = h in Cl, 


Ion 


= f a , for all |a| < to — 1 , 


I w p , — \\h\\v p , . 
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may be solved by solving the Dirichlet problem with homogeneous boundary data 

Lw = h — LF in fl, 
d a w\ gn = 0 for all \a\ < m — 1, 


\w p . — \\h\\v p . 

m,l — s — l/p —m,l — s—l/p 


\\LF \\ V , 

— m,l — s—l/p 


for some extension F and then letting u = w + F. 

Some limited results are available in the case where the coefficients c&(x) sat¬ 
isfy no smoothness assumptions whatsoever. By E3 Theorem 5.1], the Dirichlet 
problem (15.51) . in a Lipschitz domain (whose unit outward normal need not be 
in VMO), with data in appropriate spaces, has a unique solution u that satisfies 
the estimate 


(5.9) 


\w p 


,1 — s — l/p 


<C\\f\\ WA r_ i+{dn) +C\\h\\ 


v- 




provided that |p — 2| and |s — 1/2| are small enough. 

We comment on the estimate m- First, by [4j Propositon S] (listed above 
as formula (14.39D ). if u is biharmonic then the estimate (15.71) is equivalent to the 
estimate (lOll of m- Second, by (14.241) . if the coefficients aP^ are constant, one 
can draw connections between for s = 0, 1 and the nontangential maximal 
estimates of the Dirichlet or regularity problems (14.211) or (14.221) . However, as we 
pointed out earlier, this endpoint case, corresponding to the true L p -Dirichlet and 
regularity problems, has not been achieved. 


5.3. The Dirichlet problem for operators in composition form. Let us now 

discuss variable-coefficient fourth-order operators in composition form. Recall that 
this particular form arises naturally when considering the transformation of the 
bilaplacian under a pull-back from a Lipschitz domain (cf. (12.81) 1. The authors of 
the present paper have shown the wcll-posedness, for a class of such operators, of 
the Dirichlet problem with boundary data in Lr, thus establishing the first results 
concerning the L p -Dirichlet problem for variable-coefficient higher-order operators. 
Consider the Dirichlet problem 

L*(aLu) = 0 in Cl, 

u = f on dfl, 

* v ■ AXu = g on <9fl, 

t ll^(Vu)|| L 2 (on) < C||V/|| L 2 (an) + C'llfl , |lz,: 2 (dfi)- 


Here L is a second-order divergence form differential operator L = — divH(^Y)V, 
and a is a scalar-valued function. (For rough coefficients A, the exact weak definition 
of L*(aLu) = 0 is somewhat delicate, and so we refer the reader to [24].) The 
domain fl is taken to be the domain above a Lipschitz graph, that is, Q = {{x,t) ■ 
x £ R" -1 , t > (^(x)} for some function ip with £ L°°(R" _1 ). As pointed out 
above, the class of equations L*{aLu) = 0 is preserved by a change of variables, 
and so well-posedness of the Dirichlet problem (15.101) in such domains follows from 
well-posedness in upper half-spaces R™. Hence, in the remainder of this section, 
H = R". 

The appropriate ellipticity condition is then 
(5.11) A < a(X) < A, AM 2 < RerJ t A(X)r], |A(X)| < A 
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for all X £ R" and all 77 £ C”, for some constants A > A > 0. The modified 
nontangential maximal function N(X7u ), defined by 

N(Xu)(Q)= sup (/ |Vu| 2 ) ' , 

xer(Q) \is(A,dist(x,an)/ 2 ) / 

is taken from [69 ! and is fairly common in the study of variable-coefficient elliptic 
operators. 

In this case, we say that u\ ga = / and v ■ ASJu = g if 

^lim ||u( • + te) - f\\ W 2 {m) = 0, 

( lim \\v ■ AXu( ■ + te) - g\\ L 2 {dn) = 0 

where e = e n is the unit vector in the vertical direction. Notice that by the 
restriction on the domain fl, e is transverse to the boundary at all points. We 
usually refer to the vertical direction as the f-direction, and if some function depends 
only on the first n — 1 coordinates, we say that function is f-independent. 

In [24], the authors of the present paper have shown that if n > 3, and if a 
and A satisfy (15.111) and are t- independent, then for every / £ W 2 (dVL) and every 
g £ L 2 (d fi), there exists a u that satisfies (15.101) . provided that the second order 
operator L = div AS7 is good from the point of view of the second order theory. 

Without going into the details, we mention that there are certain restrictions 
on the coefficients A necessary to ensure the well-posedness even of the of the 
corresponding second-order boundary value problems; see ]28]. The key issues are 
good behavior in the direction transverse to the boundary, and symmetry. See [62] 
[69] for results for symmetric f-independent coefficients, [67ll70l 1231 [58l 159] for wcll- 
posedness results and important counterexamples for non-symmetric coefficents, 
and mm for perturbation results for t-independent coefficients. 

In particular, using the results of Emm we have established that the L 2 - 
Dirichlet problem (15.101) in the upper half-space is well-posed, provided the co¬ 
efficients a and A satisfy (15.111) and are t-independent, if in addition one of the 
following conditions holds: 

(1) The matrix A is real and symmetric, 

(2) The matrix A is constant, 

(3) The matrix A is in block form (see Section [5761) and the Schwartz kernel 
Wt(X,Y) of the operator e~ tL satisfies certain pointwise bounds, or 

(4) There is some matrix Aq, satisfying ([]]), @ or ©, that again satisfies 
(15.111) and is t-independent, such that \\A — Ao|| L oo( R „-i) is small enough 
(depending only on the constants A, A in (15.111) 1. 

The solutions to (15.101) take the following form. Inspired by formula (14.101) (taken 
from m), and a similar representation in [117] . we let 

(5.12) £h= [ F(X,Y)-f—d 2 n S*h{Y)dY 

Jn a \ Y ) 

for h defined on dfl, where 5* is the (second-order) single layer potential associated 
to L* and F is the fundamental solution associated to L. Then a(X) L(£h)(X) = 
d 2 n SJ(X) in fl (and is zero in its complement); if A* is t-independent, then 
L*(d 2 S*h) = d 2 L*(S*h) = 0. Thus u = w + £h is a solution to (15.101) . for any 
solution w to Lw = 0. The estimate ||1V(V£ h)\\ L2 ^ dQ j < ||/i|| L 2 ( 9f2 ) must then be 
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established. In the case of biharmonic functions (considered in 1 17i ). this estimate 
follows from the boundedness of the Cauchy integral; in the case of (15.10p . this is 
the most delicate part of the construction, as the operators involved are far from 
being Calderon-Zygmund kernels. Once this estimate has been established, it can 
be shown, by an argument that precisely parallels that of inzi. that there exists a 
w and h such that Lw = 0 and u = w + £h solves (15.101) . 


5.4. The fundamental solution. A set of important tools, and interesting objects 
of study in their own right, are the fundamental solutions and Green’s functions of 
differential operators in various domains. To mention some applications presented 
in this survey, recall from Sections 13.21 and 13.31 that bounds on Green’s functions 
G are closely tied to maximum principle estimates, and from Sections 13.41 and 13.51 
that the fundamental solution F is used to establish regularity of boundary points 
(that is, the Wiener criterion). See in particular Theorem 13.471 

Furthermore, fundamental solutions and Green’s functions are often crucial ele¬ 
ments of the construction of solutions to boundary-value problems. In the case of 
boundary-value problems in divergence form, the fundamental solution or Green’s 
function for the corresponding higher-order operators are often useful; see the con¬ 
structions in [426 '. 31 . 120) 1119' , 141) 193] , or in formulas (14.81) and (14.91) above. In the 
case of operators in composition form, it is often more appropriate to use the funda¬ 
mental solution for the lower-order components; see, for example, formulas (14.101) . 
(14.261) and (15.121) . or the paper 11401 . which makes extensive use of the Green’s 
function for (—A)" 1 to solve boundary-value problems for (—A) m+1 . 

We now discuss some constructions of the fundamental solution. In the case 
of the biharmonic equation, and more generally in the case of constant-coefficient 
equations, the fundamental solution may be found in a fairly straightforward fash¬ 
ion, for example, by use of the Fourier transform; see, for example, formulas (ED 
and (13.101) above, [128) 1111 , ED 11151 [511 (the relevant results of which are sum¬ 
marized as [1091 Theorem 4.2]), or [43) HD- In the case of variable-coefficient 
second-order operators, the fundamental solution has been constructed in EH EH 
[ 68 l l5p l46l l60l 122] under progressively weaker assumptions on the operators. The 
most recent of these papers, m, constructs the fundamental solution F for a 
second-order operator L under the assumption that if Lu = 0 in some ball B(X , r), 
then we have the local boundedness estimate 

(5.13) \ u {X)\<c(\f M 2 ) ' 

\ r JB(X,r) / 

for some constant C depending only on L and not on u, X or r. This assumption 
is true if L is a scalar second-order operator with real coefficients (see [112] ) but is 
not necessarily true for more general elliptic operators (see [53]). 

If Lu = 0 in B(X,r) for some elliptic operator of order 2 m, where 2m > n, 
then the local boundedness estimate (15.131) follows from the Poincare inequality, 
the Caccioppoli inequality 



(5.14) 
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and Money’s inequality 


N 

\u(X)\ <C£V 
i =o 


IB( X,r/2) 


|V J 


1/2 


whenever N/2 > n. 


Weaker versions of the Caccioppoli inequality (15.141) (that is, bounds with higher- 
order derivatives appearing on the right-hand side) were established in [3D] and 
jT71 . In [51], the first author of the present paper established the full Caccioppoli 
inequality (15.141) , thus establishing that if 2 to > n then solutions to Lu = 0 satisfy 
the estimate (15.131) . (Compare the results of Section 13.21 in which solutions to 
(—A ) m u = / are shown to be pointwise bounded only if 2 m > n — 2; as observed 
in that section, Morrey’s inequality yields one fewer degree of smoothness but was 
still adequate for the purpose of |5T].) 

Working much as in the second-order papers listed above, Barton then con¬ 
structed the fundamental solution for divergence-form differential operators L with 
order 2 m > n and with bounded coefficients satisfying the ellipticity condition (15.31) . 
In the case of operators L with 2 m < n, she then constructed an auxiliary operator 
L with 2m > n and used the fundamental solution for L to construct the funda¬ 
mental solution for L; this technique was also used in m to pass from operators of 
high order to operators of arbitrary order, and for similar reasons (i.e., to exploit 
pointwise bounds present only in the case of operators of very high order). 

This technique allowed the proof of the following theorem, the main result of 

El- 


Theorem 5.15. Let L be a divergence-form operator of order 2m, acting on func¬ 
tions defined on R n , that satisfies the ellipticity condition |53| and whose coeffi¬ 
cients A are pointwise bounded. Then there exists an array of functions F^ k (x,y) 
with the following propeHies. 

Let q and s be two integers that satisfy q + s < n and the bounds 0 < q < 
min(m,n/2), 0 < s < min(m, n/2). 

Then there is some e > 0 such that if Xq £ R", if 0 < 4r < R, if A(xq,R) = 
B(xq, 2 R) \ B(xq, R), and if q < n/2 then 

(5.16) [ [ \V n x l - s V™- q F L (x,y)f dxdy < Cr 2q R 2s 

Jy£B(xo,r) Jx£A(xq,R) 

If q = n/2 then we instead have the bound 

(5.17) [ [ \\7™- s V™- q F L {x,y)\ 2 dxdy<C{5)r 2q R 2s 

JydB(x 0 ,r) JxdA(x 0 ,R) 

for all S > 0 and some constant C(S) depending on 5. 

We also have the symmetry property 

(5.18) &].dyF^ k [x, y) = d2d 5 y Ff*. (y, x) 

as locally L 2 functions, for all multiindices 7 , 5 with |y| = m — q and |d| = m — s. 

If in addition q + s > 0, then for all p with 1 < p < 2 and p < n/{n — {q + s)), 
we have that 

(5.19) [ [ \V™- s V™- q F L (x,y)\ P dxdy<C(p)r 2n+p ( s+q - n) 

JB(xo,r) JB(xo,r) 

for all xq £ K n and all r > 0 . 
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Finally, there is some e > 0 such that if 2 — e < p < 2 + £ then V m n L extends 
to a bounded operator L p (R n ) i —> L p (M. n ). If 7 satisfies m — n/p < |^y| < m — 1 for 
some such p, then 


(5.20) 



d 2 d v F j L t k(x,y ) h k ,p(y)dy 


for a.e. x £ K” 


for all h £ L p (R n ) that are also locally in L p (M n ), for some P > n/{m — | 7 |). In 
the case of \a\ = m, we still have that 



for all h £ L 2 (R n ) whose support is not all o/ffi™. 


Here, if h £ L 


2/Tan 


), then II L h is the unique function in WT^R”) that satisfies 


£ i 

EE / = E E / d a ^a%d^n L h) k 

j=l\a\=m JQ j,k=l\u\=\0\=m JQ 


j,k =1 \ct\ = \(3\=m 1 


for all tp £ W£(R"). That is, u = H L h is the solution to Lu = div m h. The 
formulas (15.201) and (15.211) represent the statement that F L is the fundamental 
solution for L , that is, that L x F L (x,y) = 5 y (x) in some sense. 

Thus, [21j contains a construction of the fundamental solution for divergence- 
form operators L of arbitrary order, with no smoothness assumptions on the co¬ 
efficients of L or on solutions to Lu = 0 beyond boundedness, measurability and 
ellipticity. These results are new even in the second-order case, as there exist second- 
order operators L = — div AX7 whose solutions do not satisfy the local boundedness 
estimate (15. 13|) (see [97|, [58] ) and thus whose fundamental solution cannot be 
constructed as in [ 122 ] , 


5.5. Formulation of Neumann boundary data. Recall from Section lT8l that 
even defining the Neumann problem is a delicate matter. In the case of higher- 
order divergence form operators with variable coefficients, the Neumann problem 
has thus received little study. 

As discussed in Section 15.21 Agranovich has established some well-posedness re¬ 
sults for the inhomogeneous problem Lu = h with homogeneous Neumann bound¬ 
ary data. He has also provided a formulation of inhomogeneous Neumann boundary 
values; see 17] Section 5.2]. 

This formulation is as follows. Observe that if the test function ip does not have 
zero boundary data, then formula <E3D becomes 


(5.22) 


E 


j =1 


(■ Lu)j pj 


(-!) m E E [ d a vAX)a?^{X)dPu k {X)dX 

j,k =1 MH/81=m n 

m— 1 l p 

+ E E / B °m-1 -i ud lFi d<7 

i=o j =1 J on 


where BiU is an appropriate linear combination of the functions d a u where |a| = 
m + i. The expressions BiU may then be regarded as the Neumann data for u. 
Notice that if L is a fourth-order constant-coefficient scalar operator, then Bq = 
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— Ma and B\ = Ka, where I\a , Ma are given by (14.dill . Agranovich provided 
some brief discussion of the conditions needed to resolve the Neumann problem 
with this notion of inhomogeneous boundary data. Essentially the same notion of 
Neumann boundary data was used in |109] (a book considering only the case of 
constant coefficients); an explicit formula for BiU in this case may be found in [1091 
Proposition 4.3]. 

However, there are several major problems with this notion of Neumann bound¬ 
ary data. These difficulties arise from the fact that the different components BiU 
may have different degrees of smoothness. For example, in the case of the bihar¬ 
monic L p -Neumann problem of Section l4~8l the term M p u = —Bqu is taken in the 
space L p (dCl), while the term K p u = B±u is taken in the negative smoothness space 

wi^oci). 

If Cl is a Lipschitz domain, then the space IT ] 7 (dCl) of functions with one degree of 
smoothness on the boundary is meaningful, and so we may define W^^dCl), 1 /p + 
l/q= 1, as its dual space. However, higher degrees of smoothness on the boundary 
and thus more negative smoothness spaces W^ k (dCl) are not meaningful, and so 
this notion of boundary data is difficult to formulate on Lipschitz domains. (This 
difficulty may in some sense be circumvented by viewing the Neumann boundary 
data as lying in the dual space to WA^ n _ x+a {dCl)\ see [ 10811109] . However, this 
approach has some limits; for example, there is a rich theory of boundary value 
problems with boundary data in Hardy or Besov spaces which do not arise as dual 
spaces (i.e., with p < 1 ) and which is thus unavailable in this context.) 

Furthermore, observe that as we discussed in Section []J a core result needed 
to approach Neumann and regularity problems is a Rellich identity-type estimate, 
that is, an equivalence of norms of the Neumann and regularity boundary data of 
a solution; in the second-order case this may be stated as 

II^T M lli 2 (an) ~ \\ v ' ^^ u llL 2 (on) 

whenever divAVu = 0 in fl, for at least some domains Cl and classes of coeffi¬ 
cients A. In Section [5.61 we will discuss some possible approaches and preliminary 
results concerning higher-order boundary value problems, in which a higher-order 
generalization of the Rellich identity is crucial; thus, it will be highly convenient to 
have notions of regularity and Neumann boundary values that can both be reason¬ 
ably expected to lie in the space L 2 . 

Thus, it is often convenient to formulate Neumann boundary values in the fol¬ 
lowing way. Observe that, if Lu = 0 in Cl, and dCl is connected, then for all nice 
test functions tp, the quantity 

X X / d ° ^:i d' 3 Uk 

j,k =1 |a| = |/3|=m ^ 


depends only on the values of V m 1 <p on dC 2; thus, there exist functions Myf'u 
such that 


(5.23) 


X X [ a % 90U k 

j,k= 1 \a\=\p\=m Jn 


X 

|7|=m — 



M^u d 1 ipj da. 
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We may then consider the array of functions Mau to be the Neumann boundary 
values of u. This formulation only requires dealing with a single order of smooth¬ 
ness, but is somewhat less intuitive as there is no explicit formula for Mau(X) in 
terms of the derivatives of u evaluated at X. 

Also observe that, if we adopt this notion of Neumann boundary data, it is more 
natural to view the Dirichlet boundary values of u as the array {9 7 u| 0f2 : |q| = 
m — 1}, and not {<9 7 it| an : |y| < m — 1}, as was done in Sections 111 and [5721 The 
natural notion of regularity boundary values is then {X T d 1 u\ dn : |y| = m — 1}; 
again, all components conveniently may then be expected to have the same degree 
of smoothness. 

5.6. Open questions and preliminary results. The well-posedness results of 
Section Hu2l cover only a few classes of elliptic differential operators and some special 
boundary-value problems; the theory of boundary-value problems for higher-order 
divergence-form operators currently contains many open questions. 

Some efforts are underway to investigate these questions. Recall from m (see 
(15.91) above) that the Dirichlet problem 

(5.24) Lu = 0 in 12, d = f a , 1 _ a _ 1/p — ^'W-f\\wA'^ n _ 1+B (dQ) 

and the Poisson problem 

(5.25) Lu = h in Cl, d a u\ = 0, |M| w p <C\\h\\ V P 

are well-posed whenever \p — 2| and |s — 1/21 are small enough. A similar argument 
should establish validity of the Neumann problem 

(5.26) Lu = 0 in ft, M J / 1 u = g jn , \\ u \\w^ x _ B _ x/v ^ c llffll(WA^_ 1+s ( 3 n))* 
or Poisson problem with homogeneous Neumann data 

(5.27) Lu = h in 12, Mau = 0, IM|t^p < C||/i|| v p 

m,l —s —1/p -m,l-s-l/p 

Turning to a broader range of exponents p and s, we observe that perturbative 
results for the Poisson problems (15.251) and (15.271) are often fairly straightforward 
to establish. That is, with some modifications to the relevant function spaces, it is 
possible to show that if (15.251) or (15.271) is well-posed in some bounded domain 12, 
for some operator L 0 and for some l<p<oo,0<s<l, and certain technical 
assumptions are satisfied, then the same problem must also be well-posed for any 
operator L\ whose coefficients are sufficiently close to those of Lq (in the L°° norm). 

Recall from Section l5T2l that for any 1 < p < oo and 0 < s < 1, the Dirichlet prob¬ 
lem (|5.24l) for boundary data / in the fractional smoothness space WA I ^ n _ 1+s {d^l), 
can be reduced to well-posedness of the Poisson problem H5.25[l . (Some results are 
also available at the endpoint p = oo and in the case p < 1; the integer smoothness 
endpoints s = 0 and s = 1 generally must be studied using entirely different ap¬ 
proaches.) A similar argument shows that well-posedness of the Neumann problem 
(15.26P follows from well-posedness of the Poisson problem (15.271) for 1 < p < oo. (In 
the case of the Neumann problem results for p < 1 are somewhat more involved.) 
A paper (201 containing these perturbative results is currently in preparation by 
the first author of the present paper. 
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A key component in the construction of solutions of |20j are appropriate layer 
potentials, specifically, the Newton potential and the double and single layer po¬ 
tentials given by 


(n L h)j{x) = / d y F j\k( x ^y) h kAy)dy, 

k= 1 |ct|=rn R " 

£ r 

(Pnf)i(x) = iQc(x)fi(x) — Y H /_ d y F i^j( x ^y) a a/}(y) dl3 My)dy, 

|a|=|,S|=m j,fc=l 


= Y / d J F j L A x ’V) 9k^(,y) dcr(y) 

fe=l| T |=m-l , ' an 


where denotes the fundamental solution discussed in Section 15.41 and where 
/ is any function that satisfies d 1 fk = fk,~i on dQ. We remark that these are 
very natural generalizations of layer potentials in the second-order case, and also of 
various potential operators used in the theory of constant coefficient higher order 
differential equations; see in particular 1109 . 

In particular, to solve the Dirichlet or Neumann problems (15.241) or (15.261) . it was 
necessary to establish the bounds on layer potentials 


W, 


nf\\w p 1 u ^CWfWwA^ 

m,l — s — l/p rr 


ISn 


In 


l + 3 (^)’ 

these bounds are derived from the bound 


sll 


w* 


,1 — s — 1/p 


<C\\g\\ 

BP>Ps-l(dfi) 


in- 


J h\\w p , — C\\h\\ Ki ■ 

m,l — s — l/p 0,1 — s — 1/p 


This bound is the technical assumption mentioned above; we remark that it is 
stable under perturbation and is always valid if p = 2 and s = 1/2. 

Analogy with the second-order case suggests that, in order to establish well- 
posedness of the L 2 -Dirichlet, L 2 -Neumann and L 2 -regularity problems, a good 
first step would be to establish the estimates 

\\v£f\\ x < c'li/n L 2 (9n)> ll^n/Hg < c\\ ^ rf\\ L 2( dn y ll^lls) < C\\g\\ L 2 (0f2) 


for some spaces X and 2). (We remark that the corresponding bounds for constant 
coefficient operators are Theorem 4.7 and Proposition 5.2 in [109] . and therein were 
used to establish well-posedness results.) In the paper [22], now in the later stages 
of preparation, Steve Hofmann together with the authors of the present paper have 
established the bounds 


(5.28) 



\X7 m d t S A g{x,t)\ 2 tdtdx < C\\g\\ 2 L2(dn) 


(5.29) [ \V m d t V A f( X ,t)\ 2 tdtdx<C\\V x f\\ 2 L2(an) 

JR «■ Jo 

where V A = T >^ n , S A = S^ n , for scalar operators L , provided that the coeffi¬ 
cients a a p are pointwise bounded, elliptic in the sense of (15.31) . and are constant 
in the t-direction, that is, the direction transverse to the boundary of R". As dis¬ 
cussed in Section f5.3l this assumption of f-independent coefficients is very common 
in the theory of second-order differential equations. 

We hope that in a future paper we may be able to extend this result to do¬ 
mains of the form Q = {(cc,f) : t > for some Lipschitz function ip; in the 
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second-order case, this generalization may be obtained automatically via a change 
of variables, but in the higher-order case this technique is only available for equa¬ 
tions in composition form ('Section 15.31) and not in divergence form. 

In the case of second-order equations — div AS7u = 0, where the matrix A of 
coefficients is real (or self-adjoint) and t- independent, a straightforward argument 
involving Green’s theorem establishes the Rcllich identity 

|| V x Zl( • , 0)|| i 2 (g R n) ~ ||— e • ^4Vlt|| i 2 (g R n) 

where — e is the unit outward normal to R”; that is, we have an equivalence of 
norms between the regularity and Neumann boundary values of a solution u to 
div dVu = 0. Together with boundedness and certain other properties of layer po¬ 
tentials, this estimate leads to well-posedness of the L 2 -regularity and L 2 -Neumann 
problems; it is then a straightforward argument to derive well-posedness of the 
Dirichlet problem. See [ 8 ], [13], [59], [231, and others. 

We have hopes that a similar argument will yield the higher-order Rellich identity 

l|V*V m 1 u( ■ , 0)|| i 2 ( 9 R n) ~ ||A'fAU|| i 2 ( 9R ».) 

where the Neumann boundary values Mau are as in Section 15.51 for solutions u 
to divergence-form equations with t-independent and self-adjoint coefficients (that 
is, coefficents that satisfy a a n = and that a similar argument will imply 

well-posedness of higher-order L 2 boundary value problems. 

The results of Section IQ may also lead to well-posedness of L 2 boundary-value 
problems for a different class of operators, namely, operators of block type. Again, 
this argument would proceed by establishing a Rcllich-type identity. 

Let us review the theory of second-order divergence-form operators L = — div AV 
in R. n+1 , where A is an (n + 1) x (n + 1), t-independent matrix in block form; that 
is, Aj in+ i = A n+ ij = 0 for 1 < j < n, and A„ + i,n-i-i = 1 . It is fairly easy 
to see that one can formally realize the solution to Lit = 0 in M" +1 , u| Rn = /, 

as the Poisson semigroup u(x,t) = e -tv ^/(x), (x,t) £ M" +1 . Then the Kato 
estimate (El essentially provides an analogue of the Rellich identity-type estimate 
for the block operator L, that is, the /^-equivalence between normal and tangential 
derivatives of the solution on the boundary 

||<9 t u(-,0)|| 

L 2 ( R”) ~ l|V*u(-,0)|| L 2 (R n ) * 

Boundedness of layer potentials for block matrices also follows from the Kato esti¬ 
mate. 

Following the same line of reasoning, one can build a higher order “block-type” 
operator L, for which the Kato estimate (15.41) of Section I5J1 would imply a certain 
comparison between normal and tangential derivatives on the boundary 

ll^t U ( ' > 0)IIl 2 (R'>) ~ ||V£ u( ■ , 0)|| L 2( R r>). 

It remains to be seen whether these bounds lead to standard well-posedness results. 
However, we would like to emphasize that such a result would be restricted to very 
special, block-type, operators. 
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